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INTRODUCTION 



Galaxies' stellar, gas, and metal content determine the ma- 
jority of their observable properties across all wavelengths. 
Hence understanding the co-evolution of these basic con- 
stituents is at the heart of developing a comprehensive the- 
ory for the formation and evolution of galaxies. Advanc- 
ing observations have now characterised these properties in 
galaxies back into the peak epoch of cosmic star formation 
and beyond. Such observations provide stringent tests for 
galaxy formation models, and offer new opportunities to 
constrain the physical processes that govern galaxy evolu- 
tion. 

Observations have revealed tight correlations between 
stars, gas, and metals in galaxies. One example is the 
relationship between stellar mass (M,) and star for- 
mation rate (SF R), called the main sequence for star- 
forming galaxies l|Noeske et al.l |2007| ): it is slightly sub- 
linear and evolves with redshift roughly independently of 



ABSTRACT 

We use cosmological hydrodynamic simulations to investigate how inflows, star for- 
mation, and outflows govern the the gaseous and metal content of galaxies within a 
hierarchical structure formation context. In our simulations, galaxy metallicities are 
established by a balance between inflows and outflows as governed by the mass out- 
flow rate, implying that the mass-metallicity relation reflects how the outflow rate 
varies with stellar mass. Gas content, meanwhile, is set by a competition between 
inflow into and gas consumption within the interstellar medium, the latter being gov- 
erned by the star formation law, while the former is impacted by both wind recycling 
and preventive feedback. Stochastic variations in the inflow rate move galaxies off the 
equilibrium mass-metallicity and mass-gas fraction relations in a manner correlated 
with star formation rate, and the scatter is set by the timescale to re-equilibrate. The 
evolution of both relations from z = 3 — >■ is slow, as individual galaxies tend to 
evolve mostly along the relations. Gas fractions at a given stellar mass slowly decrease 
with time because the cosmic inflow rate diminishes faster than the consumption rate, 
while metallicities slowly increase as infalling gas becomes more enriched. Observa- 
tions from z ^ 3 ^ are better matched by simulations employing momentum-driven 
wind scalings rather than constant wind speeds, but all models predict too low gas 
fractions at low masses and too high metallicities at high masses. All our models re- 
produce observed second-parameter trends of the mass-metallicity relation with star 
formation rate and environment, indicating that these are a consequence of equilib- 
rium and not feedback. Overall, the analytical framework of our equilibrium scenario 
broadly captures the relevant physics establishing the galaxy gas and metal content in 
simulations, which suggests that the cycle of baryonic inflows and outflows centrally 
governs the cosmic evolution of these properties in typical star-forming galaxies. 

mass (|Elbaz et all l2007l : iDaddi et al.] l2007l : iDavel 120081 '). 

Another example is the correlation between stellar mass 
and gas-phase m etallicity, called the m ass-metallicity rela- 
tion (MZR; e.g. iTremonti et al.l l2004l l. which sh ows a re- 
markably low scatter across a wide range in mass (|Lee et al.l 
I2OO6I : IZhao, Gao, fc Gull2010l '). Moreover, departures from 
the mean relatio n are strongly correlated with other 
galaxy properties jEllison et al. l2008al: ICooper et al.l 
F-eeples. Pogge . fc Stanekl 120091 : 



200s ; 



UP >^^^^^^.. ,Lara-L6pez et al.l 120101 : 

Mannucci ct alT 2010l ). The scatter i n metallicity is tight - 
est when correlated with stellar mass (|Tremonti et al.ll2004l '). 
suggesting that stellar mass is primarily responsible for gov- 
erning the metal content of galaxies. The gas content is more 
difficult to measure because all phases must be accounted 
for (atomic, molecular, and ionised), but still shows a fa irly 
tight anti-correlation with ste llar mass (jCatinella et al.l 
l20ld : iPeeples fc Shankad l2010t ). These relations evolve 
with redshift, towards lower metallicity ([Savaglio et al.l 
I2OO5I : lErb et al.|[2006al : IZahid. Kewlev. fc Bresolinll201ll ) and 
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higher gas content (|Erb et al.ll2006bl : iTacconi et al.lboid ') at 
higher redshift. 

These trends are quahtatively consistent with the 
canonical scenario for galaxy formation in which galaxies 
start out with a gaseous halo that cools onto a central disk, 
forms stars, and self-enriches wh ile consuming its gas (e.g. 
iRees fc Ostrikej|l977l : IWhite fc R ccs 1978). However, more 
detailed observations have shown that the gas consump- 
tion rates would exhaust the gas supply quickly, and there- 
fore continual replenishm ent of gas appears to be required 
to su stain star formation (jTacconi et al.l[2010l : iGenzel et al.l 
I2OIO : Papovich ct al. 2oi3). Moreover, the metallicity evolu- 
tion is quite slow, straining self-enrichment models (e.g. the 
classical G-dwarf problem). Hierarchical structure formation 
generically predicts inflow and merging that spur galaxy 
growth. But unfettered inflow grossly overproduces global 
star forma t ion, k nown as the overc ooling problem (e.g. 
iDave et all I2OOII : iBaloeh et all l200lh . Therefore feedback 
processes must strongly regulate galaxy growth. Such feed- 
back processes are expected to manifest themselves in the 
evolution of the mass, metal, and gaseous content of galax- 
ies. Hence understanding the origin and evolution of scaling 
relations between these constituents provides key insights 
into accretion and feedback processes that govern galaxy 
growth. 

Cosmological hydrodynamic simulations have advanced 
rapidly over the past decade, to a point where they can plau- 
sibly match a wide range of properties of galaxies and the in- 
tergalactic medium (IGM) across cosmic time. One recently- 
explored physical process that greatly improves concor- 
dance with observations is strong and ubiquitous galactic 
outflows. These outflows are powered by supernovae, stel- 
lar winds, and/or photons from young stars, i.e. they re- 
sult from the star formation process itself, leading to self- 
regulated growth. Qualitatively, observations indicate that 
galaxy formation must be increasingly suppressed towards 
small masses. Outflows are now direct ly observed in most 
star-forming galax ies at z > 1 (e.g. IWeiner et al.l |2009| : 
ISteidel et al.ll201(]| V By incorporating outflows as observed 
into simulations, it is possible to yield galaxy populations 
that significantly more closely resemble those observed. 

Recently, it has been found that simulations em- 
ploying outflow scalings as expected for mom e ntum- 
driven winds (iMurrav. Quataert. fc ThompsonI l2005l : 
I Zhang fc ThompsonI I2OI0I ) are among the most 
successful at r natching a wide range of d at a on 
galaxies (e.g. iDave. Fin lator. fc Oppenheimer I I2OO6I : 
iFinlator fc Pavel l2008l: lOp pcnhcimcr ct al. 20l3) and 
the IGM (e.g. lOppenheimer fc Pavel |2006| . \2o6i . |2009| : 
lOppenheimer. Pave, fc Finlatoil booj ). although notable 
discrepancies remain. These scalings assume that the mass 
outflow rate scales inversely with galaxy circular velocity, 
providing increased suppression of star formation in smaller 
systems. Such outflows also have interesting unantici- 
pated consequences. For instance, lOppenheimer fc Pavel 
(|2008l ) found that ejected wind material more often than 
not returned into galaxies, and that this so-called wind 
recycling accretion becomes stronger at higher masses 
and d ominates the global ac cretion onto galaxies at 
z < I l|Oppenheimer et al.ll2O10l ). Furthermore, winds also 
have a "preventive" feedback effect particularly in smaller 
galaxies, by adding energy to surrounding gas which 



curtails inflow into the interst ellar medium of galaxies 
relative to inflow i n to the halo (lOppenheimer et aL 2010l: 
van de Voort et al] 1201 ll : iFaucher-Giguere. Keres. fc Mai 



2OIII ). The high outflow rates and frequent re-accretion 



suggest a continual cycling of baryons between galaxies and 
their surrounding IGM, and that this cycle is responsible 
for governing the observable properties of both. 

In this series of two papers, we investigate the 
way in which inflows and outflows within a hierarchi- 
cal structure formation context govern the main con- 
stituents o f galaxies, name ly s t ars, g as, and metals. In 
iPave. Opp enheimer. fc Finl atoil l|201ll . hereafter Paper I) 
we focused on stellar masses and star formation rates. We 
argue that many of the trends seen in simulations can 
be understood within the framework of galaxies living in 
a slowly-evolving equilibrium between inflow, outflow, and 
star formation. The inflow is at early epochs supplied pri- 
marily from the (relatively) pristine IGM, while at later 
times wind re cycling brings back gas in a mass-dependent 
fashion. As in lOppenheimer et al.l l|20ld ). we showed that 
outflows produce three-tiered stellar mass and star forma- 
tion rate functions, where the middle tier is established by 
the onset of differential (i.e. mass-dependent) wind recy- 
cling. The evolution and mass dependence of the speciflc 
star formation follows trends arising from the mass accre- 
tion rate into halos, modulated by outflows. We further 
examined the satellite galaxy population, and found that 
in models they are not particularly more common or more 
bursty than central galaxies at a given mass, and that the 
main difference versus centrals is that satellites have in- 
creasingly suppressed star formation to small masses. We 
showed that momentum-driven wind scalings provide the 
best overall flt to available observations, but the agree- 
ment is only good in the range of ~ (0.1 — 1)1/*. At lower 
masses, star formation in dwarfs seems to occur at too early 
epochs in the models, and at higher masses some additional 
mechanism is required to quench star formation in mas- 
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sive galaxies (e.g. black hole feedback; Pi Matteo et al. 
Croton et al. 2006; Pe Lucia et al. 2006; Bower et a l.| 
Fprit anot ct al. 2007; Somcrvillc ct al. 2008; Gabo r et al.1 
2OIII ). OveraU, comparing these simulations to observations 



helps constrain the way in which inflows and outflows work 
together to govern the growth of galaxies' stellar component, 
while highlighting key failures of current models. 

In this paper. Paper II in this series, we extend the 
analysis of our suite of cosmological hydrodynamic simula- 
tions with outflows to examine galaxy metallicities and gas 
fractions. The primary goal is to understand how outflows 
govern scaling relations between these quantities and their 
stellar content. We will show that the equilibrium scenario 
introduced in Paper I also provides the basic intuition for 
understanding gas and metal growth. We outline a simple 
analytic formalism that captures the main features of the 
simulation results. In it, the metallicity of galaxies is gov- 
erned primarily by outflows with a secondary effect from 
enriched infall, while the gas content is governed by a com- 
petition between cosmological gas supply and the gas con- 
sumption rate set by the star formation law. Both the metal- 
licity and gas fraction are driven by cosmic inflows, which 
diminish rapidly with cosmic time, and fluctuate on shorter 
timescales resulting in deviations from the mean relations 
that correlate with star formation. By comparing to obser- 
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vations we find that the momentum-driven wind scahngs 
provide the best match to data among the models examined 
here, but once again there are significant discrepancies at 
the highest and lowest masses. These results highlight how 
galactic outflows are a key moderator of the stellar, metal, 
and gas content of galaxies at all epochs, and in turn obser- 
vations of these properties provide valuable insights into the 
cosmic ecosystem within which galaxies form and grow. 

The paper is organized as follows. In 32] we briefly de- 
scribe our hydrodynamic simulations including our galactic 
outflow models. In ^we examine simulated mass-metallicity 
relations across cosmic time, and present a simple framework 
for understanding their physical origin. In ^ we similarly 
examine galaxy gas fractions. In ^we compare to observa- 
tions of metal and gas content to identify broad constraints 
on feedback processes. In !j6] we discuss second-parameter 
dependences of the MZR and gas fractions with star forma- 
tion rate and environment. In i}7]we explore how individual 
galaxies evolve in the mass-metallicity and mass-gas frac- 
tion planes. Finally, we summarize and discuss the broader 
implications of our work in [JH] 



2 SIMULATIONS 

The suite of simulations employed are identical to that in 
Paper I. We refer the reader there for a full discussion of all 
details, and here briefly review some of the key aspects. 



2.1 Runs 

Our simulations are run with an extended version of the 
Gadget-2 N -body -I- Smoo thed Particle Hydrodynamic 
(SPH) code (Springcl 20051). We assume a ACDM cos- 
mology (,Hinshaw ct al. 20091): = 0.28, Qa = 0.72, 
h = Ji/o/(100 kms^^ Mpc~^) — 0.7, a primordial power 
spectrum index n = 0.96, an amplitude of the mass fluc- 
tuations scaled to erg — 0.82, and Sib ~ 0.046. We call 
this the r-series, where our general naming convention is 
T[boxsize]n[particles/side][wind model]. Our primary runs use 
a boxsize of 48/i~^Mpc on a side with 384^ dark matter and 
384'^ gas particles, and a softening length of e = 2.5/i~^kpc 
(comoving, Plummer equivalent). To expand our dynamic 
range we run two additional sets of simulations with 2 x 256"^ 
particles identical to the primary runs, except one having a 
boxsize of 24/!,~^Mpc and e = 1.875/i~^kpc, and the other 
with a boxsize of 48/i~^Mpc and e = 3.75/i~^kpc. SPH parti- 
cle masses are 3.6xlO^M0, LSxIO'^Mq, and 12xlO'^M0 for 
the r48n384, r24n256, and r48n256 series, respectively, and 
dark matter particles masses are approximately 5x larger. 

Our version of Gadget-2 includes cooling pro- 
cesses us ing the primordial abundances as described in 
Katz et al . (1996)_and metal-line cooling as described in 



of cold clouds embedded in a warm ionized medium foUow- 



Oppenheimer fc Pavel ^200^ ). We inclu de heating from 



ajnet agalactia photo-ionising flux from iHaardt &: Madaul 
(|200ll ). but this has essentially no effect on galaxies that 
we can resolve since their halo masses are well above 
the mass where photo-ionisation strongly suppre sses galaxy 
formation (e.g. lOkamoto. Gao. fc Theung l2008t ). Star for- 
mation is modeled u s ing a subgrid recipe introduced in 
ISpringel fc Hernguistl l|2003l ) where a gas particle above a 
density threshold of nu = 0.13 cm~^ is modeled as a fraction 



ing lMcKee fc OstrikeJ (19771) . Star formation (SF) follows 
a Schmidt law (|Schmidtlll959l ) with the SF timescale scaled 



to m atch the z = Kennicutt law (|Kennic utt 1998). We 
use a IChabrieil (|2003l ) initial mass function (IMF) through- 
out. We account for metal enrichment from Type II su- 
pernovae (SNe), Type la SNe, and AGB stars, and we 
tr ack four elements (C,0,Si,Fe ) individually, as described 
in IOppenheime~ fc Pavel (|2008D . 

Galactic outflows are implemented using a Monte Carlo 
approach analogous to star formation. Outflows are directly 
tied to the SFR, using the relation Mwind = rjMsF, where 77 
is deflned as the mass loading factor. The probabilities for a 
gas particle to spawn a star particle are calculated from the 
subgrid model described above, and the probability to be 
launched in a wind is times that. If the particle is selected 
to be launched, it is given an additional velocity of v^^ in the 
direction of vxa, where v and a are the particle's instanta- 
neous velocity and acceleration, respectively. Choices of the 
parameters 77 and deflne the "wind model". Once a gas 
particle is launched, its hydrodynamic (not gravitational) 
forces are turned off until either 1.95 x 10^'^/(uw( kms~^)) 
years have passed or, more often, the gas particle has reached 
10% of the SF critical density. This attempts to mock up 
chimneys generated by outflows that allow relatively unfet- 
tered escape from the galactic ISM, and which are not prop- 
erly captured by the spherically-averaging SPH algorithm 
at >kpc resolution; it also yi elds results that are less sens i- 
tive to numerical resolution jSpringel fc Hernguistl [2003bh . 
For a further discussion of hy drodynamic decoupling, see 
iPalla Vecchia fc Schav^ ('2008') and Paper I. 
For this paper we run four wind models: 

(i) No winds (nw), where we do not include outflows (i.e. 

(ii) Constant winds (cw), where rj — 2 and v„ = 
680 kms~^ for all galaxies; 

(iii) Slow winds (sw), where 77 = 2 and Uw ~ 340 kms~^ 
for all galaxies; and 

(iv) Momentum-conserving winds (vzw), where galax- 
ies are identified on-tlic-fiy and their velocity dispersion a is 
estimated fsee lOppenheimer fc Pavell2008l ). and then 



(To 

V ^ —, 
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(1) 
(2) 



where /l = [1.05, 2] is the luminosity factor in units of 
the galactic Eddington luminosity (i.e. the critical lumi- 
nosity necessary to expel gas from the galax;y potential), 
and (To ~ 150 km/s is the normalization of the mass load- 
ing fact or. Choices for the former are t aken from obser- 
vations (|Rupke. Veilleux fc Sandersllioosi ). while the latter 
is broadl y is constrained to match h igh-redshift IGM en- 
richment (jOppenheimer fc Pavj[2008l ). The velocity disper- 
sion a is estimated from the bar yonic galaxy mass i\fgai us- 
ing (jOppenheimer fc Pavell2008t ) 



a = 200 



gal 



5 x IO12/1-IM0 Qi Ho 



j km/s. (3) 



See Paper I for further details. We particularly note §2.2 
for a discussion about issues related to the momentum bud- 
get in the vzw model. To reiterate, while the energy bud- 
get of this model is well within that provided by super- 
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novae, the momentum required to eject gas in the vzw 
model significantly exceeds the photon momentum input 
from young stars. A physically realistic model would then 
need to either invoke gas that is optically thick in the far 
infrared in order to extract additional momentum from pho- 
tons after reprocessing by dust (e .g. as in the models of 
iHopkins. Quataert. fc Murravl201ll 'l , or else that it is a com- 
bination of supernovae and young stars that drives the out- 
flow (e.g. as forwarded by iMurrav. Menard, fc ThompsonI 
I2OIOI ). 

2.2 Computing galaxy metallicities and gas 
fractions 

We use SKIE0 (Spline Kernel Interpolative Denmax) to 
identif y galaxies as boun d groups of star-forming g as and 
stars ^Keres et all l2005l : lOppenheimer et al.l I2OI0I ). Our 
gala xy stellar mass lim it is set to be ^ 64 star parti- 
cles (|Finlator et al.|[2006l ). resulting in a minimum resolved 
mass of Mgal = 1.1 X 10^ Mq in our r48n384 series of runs. 
We will only consider galaxies with stellar mass M, ^ 
Mgai in our analysis. We separate galaxies into central and 
satellite galaxies by associating each galaxy with a halo, 
where we identify halos via a spherical overdensity algo- 
rithm ijKeres et al.ll2005l ). 

To compute gas fractions, we set Mgas to be the mass 
of all star-forming gas in a SKID-identified galaxy. We then 
define 



where M* is the stellar mass of the galaxy. Note that 
some authors choose /ga s ~ M^i^^/M,, and for instance, 
iPeeples fc Shankad (|201( jl argue that this definition is more 
natural in terms of understanding the origin of the MZR. 
Nevertheless, we prefer the above definition because it intu- 
itively translates into the fraction of a galaxy's (baryonic) 
mass that is in gas. In the end, so long as models and data 
are compared using the same quantity, the exact definition 
is not critical. 

Our gas mass includes all phases of the ISM, including 
the cold neutral medium, molecular gas, and the warm ion- 
ized medium. The latter typically makes a small mass contri- 
bution, but the relation between the first two depends on the 
internal physics of the ISM (particularly self- shielding) that 
our s imulations do not accurately track (see lPopping et ali 
I2OO9I . for further discussion). It is therefore important to 
compare to data that accounts for both neutral (Hi) and 
molecular {H2) components, as well as having been cor- 
rected for Helium. Furthermore, gas mass measurements can 
be sensitive to surface brightness effects in the outer regions 
of galaxies. 

The determination of gas content in simulated galax- 
ies depends on the density threshold for star formation. 
Here we assume this density threshold is nn = 0.13 cm~"^, 
which is a somewhat arbitrary choice motivated b y res- 
olution considerations (Springcl fc Hernguistl l2003bl ). This 
assumption makes little difference for galaxy star forma- 
tion histories, because dynamical perturbations are typi- 

^ http://www-hpcc.astro.washington.edu/tools/skid.html 



cally so frequent that gas is rapidly driven inwards un- 
til it forms stars, and thus raising (lowering) the thresh- 
old would only introduce a small dela y (ad vancement) in 
processing gas into stars (jSchave et alTlioiol ). Naively, gas 
fractions would be more affected since changing the thresh- 
old density will include more or less gas as "star-forming". 
But it is not straightforward to even determine the sign of 
the effect. For instance, lowering the threshold would pro- 
vide more gas eligible for star formation, but would also in- 
crease the star formation rate, thus lowering the gas fraction. 
In high-resolution simulations of individual galaxies (e.g. 
iGuedes et al]|201ll ). it is seen that raising the threshold den- 
sity does lower the gas fraction, but this is also sensitively 
dependent on their feedback algorithm which is quite dif- 
ferent than ours. Hence gas fraction comparisons should be 
taken with some caution, and are mainly intended to illus- 
trate trends. 

The galaxy gas-phase metallicity is defined as the SF- 
weighted metallicity of all gas particles in SKID-identified 
galaxies. This definition most closely mimics how metallic- 
ities are measured observationally using nebular emission 
lines emanating from star-forming regions. We use the oxy- 
gen abundance as a metallicity tracer in our models, since 
this is also typical of observational determi nations. We as- 
sume a solar oxygen mass fraction of 0.00574 l|Asplund et al.l 
|2009|) . or [O/H]0 -I- 12 = 8.69. The weighting of metallicity 
by star formation mitigates the issues regarding the outer 
regions of galaxies that plague gas fractions, since star for- 
mation is typically concentrated in the central region of the 
galaxy. 



3 GALAXY METAL CONTENT 

In this section we will examine the drivers behind the mass- 
metallicity relation (MZR) and its evolution out to high red- 
shifts. We will focus on the physical mechanisms that con- 
nect inflows and outflows to the observable metal content 
of galaxies, in particular placing the form and evolution of 
the MZR within the context of the equilibrium scenario for 
galaxy evolution. 

3.1 The z — mass-metallicity relation 

Figure [T] shows the z = relation between stellar mass 
and gas-phase metallicity, the MZR, in our simulations. The 
simulation data points are color-coded by star formation 
rate within bins of stellar mass, which we will discuss in 
^ The main body of points comes from the r48n384 runs; 
the smaller-volume r24n256 run galaxies are shown as the 
sparser set of points at M, < 1.1 x 10^ to extend the dy- 
namic range. The magenta curve shows a running median, 
and the errorbars show la deviations about the median. For 
comparison, the SDSS mass-metallicity relation mean (thick 
line) and la scatter (dashed lines) is overlaid, but we will 
defer discussion of comparisons to observations until !j5l 

Metallicity measures have an uncertain normalization. 
This com es from uncertainties in metallicity determina- 
tions (e.g. lKewlev fc Ellison|[2"008l ) . u ncertainties in adopted 
metal yields (see discussion in lOppenheimer fc Pavel 
|2008| ). and uncertainties in the solar metal abun- 
dance lAsplund et all |2009| ). Hence we treat the overall 
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Figure 1. The M« — Zgas relation (MZR) at z = in our r48n384 cosmological hydrodynamic simulations employing our four galactic 
outflow scalings: momentum-driven scalings (upper left), constant winds (upper right), no winds (lower left), and slow winds (lower right). 
Coloured points represent individual simulated galaxies, colour-coded by SFR within bins of M* into upper (blue), middle (green), and 
lower (red) thirds. Magenta li nes show a running m edian of the simulated points, with la scatter about the median. The thick solid line 
is the z ^ MZR from SDSS jTremonti et al.ll2004l ) with dashed lines showing the range enclosing 16% — 84% of the data. Note that all 
model oxygen abundances have been multiplied by 0.8 in order to match the amplitude of the observed MZR at M« ^ IO^^Mq, which 
is within systematic uncertainties in metallicity measures; the shape, scatter, and evolution are independent predictions. 



metallicity normalization as a free parameter. Given that 
we do not have any form of feedback that que nches massive 
galaxies in these runs (e.g. iGabor et al.ll201ll ). our simula- 
tions most robustly model star-forming galaxies at masses 
below M* . Therefore we normalize our metallicities to the 
observed MZR at M« ^ IO^^Mq, where it so happens that 
all our wind simulations predict roughly similar metallici- 
ties. This normalization requires us to multiply all simulated 
metallicities by an arbitrary factor of 0.8. We apply this 
same factor at all redshifts and for all models. Hence the 
independent predictions of our simulations are the shape, 
slope, scatter, and evolution of the MZR, but not its overall 
amplitude. 

All the wind models produce a general trend of increas- 



ing metallicity with mass and a turnover to flat at high 
masses. The no- wind case in contrast produces a nearly 
flat MZR, as one would expect from e.g. closed box evo- 
lution. All the wind models approach the no- wind case 
at Af» > 10^^ Mq, as these wind models all eject pro- 
portionally less material from the most massive galaxies, 
an d the material that is eje cted tends to quickly fall back 
in (jOppenheimer et al.ll2010l ). If some form of ejective feed- 
back to quench massive g alaxies was included in our mod- 
els (e.g. iGabor et al]|201ll ). this plateau metallicity may be 
lower. 

While broadly similar, there are clear differences be- 
tween various wind models. The constant wind model yields 
another turnover at low masses (M, < 10^° Mq) towards a 
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flat MZR. The slow wind model produces a similar turnover 
at somewhat smaller masses. The momentum-driven scal- 
ings case does not produce such a flattening, at least within 
the mass range probed by these simulations; the MZR slope 
here is nearly constant at Z oc Af»'^ for M, < 10^°-^ Mq. 

To understand the origin of these features for 
the various win d mod els, we review the findings from 
iFinlator fc Davj ([2008') who developed a simple analytic 
understanding of the MZR. In their model, the gas-phase 
metallicity of a galaxy is set by a balance between inflow and 
outflow plus star formation. Inflow brings in low- metallicity 
gas, star formation enriches that gas, while outflows modu- 
late the fraction of inflow that turns into stars. In equilib- 
rium, the three terms are related by 

Minflow = M, + Moutflow (5) 

(see also lDutton et al.ll2O10l : iBouche et al.ll2010l }. Rewriting 
this in terms of the mass loading factor 77 = Moutflow /A/,, 
we obtain 

Minflow = (l + r?)M. (6) 

In this simple "equilibrium" picture, the metallicity is given 
by the amount of metals produced by star formation, which 
is j/M* where y is the yield of metals per unit star formation, 
divided by the rate of gas inflow to be enriched. Hence 

Z = yAf,/M„flow = (7) 

Here we have, for simplicity, assumed that the infalling gas 
has negligible metallicity; we will relax this assumption in 
37] (see Equation [TT|) . 

A key assumption in this formalism is that the mass 
loading factor reflects the amount of material that is ejected 
from the galaxy without returning quickly. In this sense, rj 
shoul d be regarded as an "eff ective" mass loading factor, 
which IFinlator fc Davj l|2008l ) showed generally tracks the 
input value of ri in our simulations when the wind velocity is 
comparable to or exceeding the escape velocity. We reiterate 
that our simulations hydrodynamically decouple outflowing 
gas; simulations that choose not do so can have significantly 
different effective r; despite h aving the same input 77 (e.g. 
iDalla Vecchia fc Schavell2008l ). The similarity between our 
input ri and effective rj (at least above the escape velocity) 
thus reflects this particular modeling choice. 

It is worth noting that Equation[7]does not have any ex- 
plicit dependence on stellar mass, but only depends on inflow 
and outflow rates. The physical interpretation is that the 
gas-phase metallicity does not reflect a historical record of 
star formation in a galaxy (as in a closed-box scenario), but 
rather reflects its recent (i.e. over a gas depletion timescale; 
>i4.2[) balance between inflows and outflows. This then distin- 
guishes stellar metallicities, which must reflect the history 
of metal buildup, from gas-phase ones. In practice, how- 
ever, the fairly slow evolution of the MZR ( H3.2p means that 
galaxy stellar metallicities are only slightly lower than gas- 
phase ones. We leave a detailed examination of stellar vs. 
gas-phase metallicities for future work. 

Using Equation [71 we can understand the behavior of 
the various wind models. In the no-wind case, 77 = 0, and 
the metallicity is therefore close to constant. Although the 
MZR is close to flat, there remains a slightly slope owing to 
the rapidity of infall and the lack of reduction of fresh gas 



by outflows, which results in galaxies being not quite able to 
attain equilibrium. We will see in the next section that this 
effect becomes exacerbated at higher redshifts. Furthermore, 
there is more enriched infall into higher mass galaxies, as we 
discuss in ijTl 

The momentum-driven wind scalings assume rj oc v^^ oc 

M* (approximately). Hence when rj ^ 1, this approxi- 

1/3 

mates Z oc A/« . The turnover at high masses is set by the 
normalization of rj, namely ao, which produces rj < 1 for 
Af, > 10" A/0 (at z = 0; this mass evolves mildly upwards 
with redshift). 

The constant and slow wind cases introduce another 
consideration: The competition between wind speed and 
escape velocity. Unlike in the momentum-driven scalings 
where v„ ~ Vcsc, here there is a transition mass above 
which the winds cannot escape, and fall quickly back into 
the galaxy. Hence the wind recycling time is very short, 
meaning that winds have little effect ([Oppenheimer et al.l 
I2OIOI ). Stated in terms of the above formalism, the effective 
mass loading factor in t hese models approach zero above 
that threshold mass (see IFinlator fc Dave'2008' and §2.2 of 
Paper I), so that at low masses ~ 2 while at high masses 
77 ~ 0, with a steep transition between these regimes across 
which the metallicity changes by a factor of 1 -I- 77 = 3. Be- 
cause the wind speed is twice as fast in the cw case, the 
transition occurs at a higher mass than in sw. Although the 
transition to the low-mass regime is not well probed at the 
dynamic range of these simulations (particularly in the sw 
case), it is still evident. This transition is also evident in 
the stellar mass and star formation rate functions in Pa- 
per I. This equilibrium scenario strongly predicts that the 
low-mass MZR will continue to be flat to small masses in 
this model. 

These simulations and the associated equilibrium model 
indicate that the MZR is critically governed by the mass out- 
flow rate from galaxies (i.e. 77) and its scaling with M«. This 
differs fundamentally from the canonical explanation for the 
MZR that invokes a competition between the galaxy poten- 
tial well and the outflow ve locity to mo dulate the me tals re- 
tained within a galaxy (e.g.lDekel fc Silk 1986; D ckel fc Wool 
I2OO3I : Irremonti et al.ll20M ): here, there is no mention of the 
potential well depth except indirectly via its effects on the 
mass loading factor. It is often canonically stated that low 
mass galaxies preferentially eject more of their metals, and 
hence have lower metallicity. In our scenario, it is not the 
ejection of metals that is modulated by outflows, it is the 
production of metals from a given amount of hierarchical 
inflow that is regulated by outflows. 

The turnover in the MZR at high masses has typically 
been thought to reflect the transition to a regime where out- 
flows can no longer escape the galaxy ([Pekel fc Silk|[l986l : 
iTremonti et al.|[2004l '). In our constant-77 cases, this is accu- 
rate. But in our momentum-driven wind scalings case, since 
the wind speed scales with the escape velocity. Hence the 
turnover is instead caused by 77 becoming less than unity at 
large M, in Equation [71 Since this model seems to predict 
an MZR that is overall in better agreement with data, par- 
ticularly to low masses, this suggests that the conventional 
interpretation of the MZR being governed by the potential 
well depth may not be accurate. 

In summary, the MZR is governed by the "star forma- 
tion efficiency" , where here we mean this in the cosmological 
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sense as the amount of infalling material that is converted 
in to starfl Thi s inter pretation agrees with the simulations 
of lBrooks et"all l|2007l ) , who also found that the MZR is gov- 
erned by the cosmological star formation efficiency that is 
modulated by supernova feedback. In our models, this star 
formation efficiency and hence the metallicity is directly con- 
trolled by the rate at which mass is ejected in outflows. 



4 GALAXY GAS CONTENT 

In this section we examine galaxy gas fractions and their 
evolution across cosmic time. As with the MZR, we attempt 
to provide physical intuition for what establishes a galaxy's 
gas fraction and its evolution by connecting it to gas inflow 
and outflow processes. 



3.2 MZR evolution 

Figure O depicts the evolution of the MZR from 2 = 3^-0 
in our four wind models. We show the running median at 
each redshift in the large panel, and the small panel below 
shows the la variance about the median within each mass 
bin. The cyan points in the lower pan els show the observed 
la variance at z ~ from SDSS data (|Tremonti et al.ll2004l ) 
for comparison; we will discuss this further in [JS] We sepa- 
rately show the results for the r48n384 (solid lines), r24n256 
(dotted), and r48n256 (dashed) runs; the good agreement 
at all overlapping masses indicates that these results are 
resolution-converged at least over the range of resolutions 
and volumes probed here. Although we don't show it here, 
all models at all epochs retain the second-parameter trend 
shown in Figure [1] wherein lower SFR galaxies at a given M* 
have higher metallicities. 

All the wind models have the general shape of their 
MZR as expected from the equilibrium model and described 
in the previous section. The shapes remain similar at all red- 
shifts, because the form of ri{M,) for each model does not 
change, and this governs the MZR shape as described in i)3.1l 
The no- wind case appears further out of equilibrium at ear- 
lier epochs, as it deviates more strongly from the expected 
behavior of a constant metallicity at all masses. This is ex- 
pected because accretion rates are higher at early epochs 
in comparison to gas processing rates within the ISM; we 
will discuss this further in i|4.3l The trend of a metallicity 
increasing with time at a given Af* is quite generic, at least 
among the models considered here. This is not trivial; it is 
possible to design models that are qui te reasonable in many 
ways but yield the opposite evolution (|Arrigoni et al. Iboiol . 
and M. Arrigoni, priv. comm.). 

In detail, the evolutionary rate at a given M* varies 
somewhat with wind model. The momentum-driven wind 
scalings produce little early evolution, and more evolution 
from 2 ~ 2 — >■ 0. The constant-w^, models have less late evo- 
lution, particularly in the slow wind case, and more rapid 
evolution at early epochs. In Paper I we saw that such trends 
are also seen in the evolution of the galaxy star formation 
rate functions (Figure 2 of Paper I), where the constant-«u, 
cases evolved less out to 2 ~ 2 compared to the momentum- 
driven scalings case. This qualitative similarity in evolution 
is consistent with the interpretation that the MZR is gov- 
erned primarily by galaxies' star formation rates, as sug- 
gested by Equation [T] We will discuss the evolution of the 
MZR further in ijT] when we study how individual galaxies 
evolve in mass-metallicity space. 



^ This is notably different than the definition of star formation 
efliiciency in the interstellar medium, which describes how quickly 
molecular gas is processed into stars. 



4.1 Gas fractions at 2: = 

Figure [3] shows the z — gas fractions in our suite of simula- 
tions as a function of stellar mass (the mass-gas relation, or 
MGR). The points are color-coded by star formation within 
a given mass bin as in Figure [1] this will be discussed in SjH] 
For comparison, the data points with errors show observed 
gas fractio ns (Hl-|-H2), correc t ed fo r Helium, compiled and 
binned bv lPeeples fc ShankaJ (|2010l ). 

At the massive end, all models show decreasing gas 
fractions with stellar mass. At lower masses, however, all 
wind models eventually deviate from this trend, display- 
ing a maximum typical gas fraction below which the gas 
fraction becomes lower to smaller masses. The no-wind case 
shows no such maximum. The mass at which this maxi- 
mum occurs appears to be related to wind recycling, i.e. 
the return of previo u sly-ej ected material back into a galaxy. 
lOppenheimer et al.l (|201Cl) showed that the recycling time 
becomes long at smaller masses, eventually exceeding a Hub- 
ble time. The mass at which the recycling time becomes 
longer than the Hubble time in each wind model is, to a 
good approximation, the mass at which the maximum gas 
fraction is seen. This suggests that the lower gas fractions at 
small masses occurs at least in part because ejected winds 
never return to these galaj cies. 

As em phasize d by van de Voo rt et~all l|201ll ') and 
lFaucher-"Giguere. Keres. fc Ma (2011), feedback adds energy 
to surrounding gas and prevents material from entering into 
smaller galaxies. Since much of the m aterial entering into 
galax ies' ISM at 2 = is recycled winds l|Oppenheimer et al.l 
|2010| ). the majority of the effect is that small galaxies do 
not re-accrete their winds. The constant wind model shows 
the highest turnover mass in /gas, while the slow wind and 
momentum-driven scalings cases occur at lower masses. The 
latter wind model also shows a somewhat slower drop-off in 
/gas to lower masses, reflecting its less steep dependence of 
recycling time with mass. 

Another possible issue is that perhaps many low-mass 
galaxies are satellites that have lower gas content owing to 
stripping or strangulation processes. However, we will show 
in i]6.3l that the satellite fraction does not increase appre- 
ciably to small masses, and we will demonstrate that the 
turn-down in gas fractions to small masses is present in both 
satellites and centrals. Hence this particular trend does not 
arise from satellites, although other interesting trends do 
that we will explore in H6.3I 

The non-monotonic behavior of gas fractions in wind 
models reflects a similar behavior, with similar characteristic 
scales, as the speciflc star formation rate (sSFR=SFR/A/*) 
examined in Paper I. There we saw that low-mass galax- 
ies had depressed sSFRs relative to an extrapolation from 
higher masses in all wind models; the no-wind case showed 
no such deviation. This trend arises because of a combi- 
nation of wind recycling, which brings extra accretion at 
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Figure 2. Large panels show the evolution of the mass-metallicity relation in our four wind models at ^ = 0,1,2,3. Lines show the 
running median within mass bins at each redshift. Solid lines show the results from our r48n384 simulations, dotted lines show r24n256 
runs, and dashed lines show r48n256 runs; the general consistency between the three indicates that the results are numerically converged, 
though the lowest resolution runs show noticeably lower metallicities. All models show an upwards evolution of metallicity at a given 
mass (at a rate that depends on wind model), although the characteristic MZR shape unique to each model does not change with redshift. 
Smaller panels below each large panel show t he ler scatter about th e median relation for all the models. Cyan points show the Icr 
scatter in observations of the z = MZR from FiVemonti et al ] l|2004l 'l. 



high masses (jOppenheimer et al.ll201oh . and preventive feed- 
back which suppresses accretion into galaxies at the lowest 
masses (e.g. Ivan de Voort et al.l l201ll l . The suppression of 
gas fractions is seen to be directly proportional to the sup- 
pression in sSFRs, which we will explain in Hence the 



discrepancies of models relative to observed sSFRs of dwarf 
galaxies noted in Paper I is directly traceable to lowered gas 
fractions in dwarfs predicted in models. 

In summary, gas fractions fall with mass at the high- 
est masses but show a turnover at low masses in all the 
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Figure 3. The relation between /gas and M» (MGR) at z = in our four winds models, with magenta lines showing the median and 
la scatter as in Figure [T] The points are color-coded by SFR within stellar mass bins: blue for upper third , green for middle third, r ed 
for bottom third. Data points show mean values as a function of M* from a compilation of observations bv lPeeples &: Shankad 1 I2OICII ). 



wind models. This turnover is not seen in observations, as we 
will discuss in ij5] To understand the origins of these trends, 
and also gas fraction evolution, let us examine an instruc- 
tive quantity for understanding gas processing in galaxies, 
namely the depletion time. 



4.2 Depletion time 

Figure U shows the gas depletion time tdep as a function 
of stellar mass in our four wind models, from z = 3 ^ 0. 
We define tdop = Mgas/SFR, which is the time that galaxy 
would take to consume its current gas supply at its cur- 
rent star formation rate. Our three sets of simulations for 
each wind model are shown to demonstrate good resolution 
convergence over the dynamic range probed. To illustrate a 
useful trend, we divide tdcp by the Hubble time tn at each 
redshift. 



Several general features of tdcp(A^*,z) are evident in 
Figure [4] First, higher-mass galaxies have shorter depletion 
times. Second, depletion times are, to first order, indepen- 
dent of wind model. Finally, tdep/^if is essentially invariant 
over most of cosmic time. The insensitivity to wind model 
and invariance when scaled to the Hubble time provide key 
clues into the physics governing tdcp, and in turn, galaxy gas 
contents. 

We can broadly understand these trends using a 
straightforward argument based on our star formation law. 
The depletion time measures how gas, once within a galaxy, 
gets consumed into stars. Our simulations model conver- 
sion of gas to stars by assuming a Kennicutt-Schmidt Law, 
which equivalently follows the relation that the star forma- 
tion rate is the gas mass divided by the dynamical time tdyn 
(of the star- forming disk), times some overall efficiency fac- 
tor that is measured to be around 2% b oth locally and in 
distant galaxies (e.g. iGenzel et al1l2010l ). Hence the deple- 
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Figure 4. Depletion time t^ep = Mgas/SFR in units of the Hubble time tn in our four wind models at 2 = 0, 1,2,3. Solid lines show 
results for the r48n384 series, dotted lines show r24n256, and dashed lines show r48n256. idcp drops to larger masses roughly as AI~^'^ 
and scales with tn, essentially independent of winds, as derived in Equation [S] 



tion time sh ould scale as the dynamic al time. In a canonical 
disk model (|Mo. Mao, fc Whitd [19981 ). the dynamical time 
evolves as tH, and hence we expect tdcp/iff to be approxi- 
mately constant, which is generally confirmed in Figure |31 
although with deviations at low masses a.t z — that we 
discuss below. 

The mass dependence of t^cp depends on the conver- 
sion rate of gas into stars. This is set by the details of inter- 
nal structure of galaxies within our simulation along with 
the assumed Kennicutt-Schmidt star formation law, which 
is SsF oc Sja'l where Ssf and Egas are the surface densities 



of star formation and star-forming gas, respectively. Hence 

tdop = Egas/SsFR OC Sgas' . 

We then employ an empirical relation measured in our 
simulations of Egas oc M*''* that we show in Figure [S] 
It does not evolve appreciably with redshift. While we do 
not explicitly show it, this relation holds reasonably well 
for all the wind models. Observations of the stellar mass 
density profile in late-type S PSS galaxies indicate E* oc 
M°-^^ (jKauffmann et alJl20o3 ). which is slightly shallower. 
Gas profiles are more difficult to measure, but at least in 
our models, star-forming gas generally traces stars. Using 
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Figure 5. Gas surface density Sgas (comoving) versus stellar 
mass M* in galaxies from our r48n384vzw simulation, at 2 = 
(blue) and 2 = 3 (red). The gas surface density is taken to be 
the gas mass divided by nR^, where R is the stellar half- mass 
radius of the galaxy. This relation is reasonably well described by 
Sgas oc M^^'^, as shown by the arrow. 



galaxies begin with a large reservoir of gas and consume 
them rapidly (e.g. Eggen. Lvnden-Bell. fc Sandag^ Il962l : 



iMaraston et al.ll2010l ). 



At 2 = 0, we see a systematic departure from the mean 
trend towards lower tdcp at low masses. Part of this owes to 
satellite galaxies that are increasingly starved of gas to low 
masses at low redshifts (Figure 6 of Paper I) , which occurs 
even in the no-wind case. In the wind models, an additional 
role is played by preventive feedback and the lack of wind 
recycling at small A/». As usual for effects involving wind 
recycling, the departure from the canonical relation occurs 
at a higher mass in cw relative to sw, and is more gradual in 
the vzw case. These effects are very mild at earlier epochs 
when wind recycling is not as common. 

The depletion time has sometimes been interpreted as 
a measure of star formation efficiency variation with galaxy 
mass, such that more massive galaxies more efficiently con- 
vert gas into stars. But in our models, the star formation 
efficiency is an input constant t hat is calibrated to match ob- 
servations of present-day dis ks jSpringel fc Hernauistl2003bl : 
lOppenheimer fc Dave|[200i) . and does not vary with galaxy 
mEiss. The trend of tdep with M* arises from the assumed 
star formation formation law, which is why it is very weakly 
dependent on details of feedback. 



this shallower slope would not significantly change our re- 
sults, but we prefer to use the simulated slope since we are 
trying to develop an analytic understanding of the simu- 
lation results. Our r24n256 and r48n256 simulations have 
slightly different amplitudes for this relation but the slope 
is identical, suggesting that the trend with A/« is insensitive 
to resolution at least over the (admittedly narrow) range 
probed by these simulations. 

Putting this together, we obtain 

idcp oc ^i^A^-°■^ (8) 

This scaling provides a good match to tdcp(Af,) in our sim- 
ulations at most epochs and masses as shown in Figure |31 
Since we assume the same star formation law in all our wind 
models, there is little sensitivity to outflows in this relation. 

We note that t^ep oc tn does not necessarily have to be 
the case - for instance, in mergers, tdyn is quite small within 
the dense central star-forming region during the starburst 
phase, meaning that tdop of such starburst galaxies should 
lie below the mean relation. The relatively tight relation of 
tdcp(Aft) in our models indicates that, in analogy with the 
tight main sequence (jpavc 2008, Paper I), mergers are not a 
dominant population at any redshift. Hence even though we 
have not explicitly used the equilibrium condition in deriv- 
ing tdcp, the connection between global virialization and the 
properties of the star-forming region in galaxies only holds 
when galaxies are in a steady-state situation. 

tdep is to first order insensitive to outflows; all mod- 
els, regardless of winds, show a roughly similar tdcp(Af,). 
This is expected since the derivation of Equation [8] does 
not involve any aspect that depends on outflows, effec- 
tively employing only the star formation law and virial 
arguments. This supports the idea of Paper I that star 
formation is supply-regulated, that is, star formation oc- 
curs in proportion to the gas available to form stars. 
This differs relative to expectations from models in which 



4.3 Gas fraction evolution 

Figure |5] shows the evolution of the MGR from z — 3 ^ 
in our simulation suite. Following Figure [2] the solid lines 
show medians with la spread for the r48n384 runs, while 
the dotted and dashed lines show likewise for the r24n256 
and r48n256 runs, respectively. Once again, good resolu- 
tion convergence is seen, as all the key features are repro- 
duced at each resolution, although we emphasize that the 
dynamic range probed here is only a factor of eight in mass. 
The smaller panels show the running la variance about the 
median, and cyan points show the observed scatter from 
iPeeples fc Shankaj l|2010l ). 

All models display a slowly falling gas fraction with 
time at a given mass, while the variances (small lower 
panels) become slightly higher with time. Hence higher- 
z galaxies are more gas-rich, in qualitative agreement 
with observations. The fundamental physics governing this 
is a competition between gas inflow and gas consump- 
tion rates. In the cold accretion paradigm, the amount 
of gas inflowing into the star forming region is propor- 
tional to the the gas inflowing at the halo virial radius, 
since cold streams ch annel material to the center of the 
halo (jPekel et al.ll2009l ). In detail, preventive feedback mech- 
anisms can reduce ISM inflow (Ivan de Voort et~aLl I2OIII : 
iFaucher-Giguere. Keres. fc Mai I2OIII ) particularly at low 
masses. 

The amount of gas entering the virial radius can be es- 
timated by the total mass accreti on rate times the baryon 
fraction, which scales as (l+z)^'^'^ (jPekel et al.ll2009l ). Mean- 
while, the gas consumption rate is given by how fast gas 
can be processed into stars or an outflow, which is given 
by tdcp/(l + v)- Since tdep oc tn, at any given mass (which 
approximately corresponds to a given rj), the consumption 
rate is given by tj^^. Since tj^^ evolves with z slower than 
(1 + z)^''^^ (e.g. tJi^ oc (1 -I- z)^'^ in the matter-dominated 
era), it is straightforward to see that the gas supply rate 
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Figure 6. Like Figure[2l only for gas fraction, showing the evolution of the MGR in our four wind models at 2 = 0, 1, 2, 3. Cyan points 
show the observed scatter from Peeples &: Shankar. 12010.) . The arrow in the no- wind panel shows a slope of —0.5, as is roughly expected 
at the high-mass end; see text. 



drops faster with time than the gas consumption rate. This 
explains why galaxies start out gas-rich but then the gas 
fraction drops as the gas consumption rate catches up. 

While the naive notion that galaxies simply "consume 
their gas" is qualitatively in accord with obse rvations, the 
rapid consumption times l|Tacconi et al.l |2010| ) imply that 
gas must be continually supplied. In our simulations, this 
indeed happens, but at a rate that becomes slower with time 
(relative to the consumption rate), causing a gas fraction 



that slowly drops. The no-wind case in Figure [S] shows a 
self-similar (in M« ) downwards evolution in /gas arising from 
this scenario. 

We can quantify these scalings and gain more insight 
into gas fractions by first using the definition of tdop to 
rewrite 



1 + isp/tdcp ' 
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where tsF = M,/SFR=l/sSFR. Thus the dependence of 
/gas on mass and redshift reflects the dependence of the ratio 
tsF/tdep on these quantities. 

Let us consider the redshift evolution first. In Paper I 
we showed that the observed tsF at A/* — W^^Mq was con- 
sistent with following the trend predicted by halo accretion, 
namely tsF oc (1 + z)~^''^^, from z ~ 2 ^ 0. Figure |4] shows 
that tdep oc tH over that time, albeit with some deviations 
at low-z at low masses. Hence the ratio tsp/idcp rises with 
time, approximately as (1 + z)"'^^ at high-z and (1 + z) at 
low-z, resulting in a dropping gas fraction. Since tsF and 
tdep are to rough order independent of winds, this explains 
the basic behavior of /gas dropping with time in all models. 

Equation|9]can also be used to gain insights on the mass 
dependence of /gas- At larger masses (e.g. M, > 10'°Mq), 



tdep <C tsF, so we can approximate /g. 



idepsSFR. Con- 



sider first the no-wind case, here sSFRcx M^T (approx- 
imately; see Figure 2 of Paper I). Combined with tdcp oc 
M^"'^ (Equation |S}, this then roughly predicts /gas oc 
^-0.5 'pjjjg giope is indicated in the nw panel of Figure |6l 
and provides a good fit to the high-mass slope of /gas(Af«). 
At lower masses, once t^cp/tsF becomes a significant fraction 
of unity, the gas fraction levels off. 

In the wind models, /gas will reflect features seen in the 
sSFR, and at z = there are additional deviations owing to 
tdep. In general, sSFRs in the wind models show a turnover 
to lower sSFRs at low masses (Figure 3 of Paper I). The loca- 
tion and strength of that turnover depends on the particular 
wind model, owing to wind recycling as discussed in Paper I. 
At high masses, sSFR in wind models is raised over the no- 
wind case by wind recycling, which is is more rapid in higher 
mass galaxies. At low masses, the sSFR is lowered owing 
to preventive effects of winds adding e nergy to surround- 
ing g as (e.g. lOppenheimer et al] I2OI0I : Ivan de Voort et al] 
[20l3). These trends are directly reflected in /gas(Af«): At 
high masses, they are somewhat larger than the no-wind 
case, while they all show a turnover to low gas fractions 
at low masses. At z = 0, this reduction at low masses is 
exacerbated by the drop in t^cp to low masses. 

An interesting regime that is not probed here is the 
very high-z epoch (z > 3). At sufficiently high redshifts, 
the rapidly-rising accretion rate will begin to exceed the less 
rapidly-rising gas consumption rate. In that case, star for- 
mation cannot keep up with the gas supply, and the galaxies 
will no longer be in equilibrium, this is then the gas accu- 
mulation phase. The exact epoch where this happens de- 
pends on feedback; when outfiows are highly mass loaded, 
the amount of infalling gas that needs to be processed into 
stars is reduced, and equilibrium can occur earlier on. Also, 
there is some mass dependence because idep has a significant 
mass dependence (Figure |4|, b ut tsF (or sSFR) is ess entially 
independent of mass at high-z ([Gonzalez et al.ll2010l ). Inter- 
estingly, there are now empirical constra ints on the gas ac- 
cumulation epoch: IPapovich et al.l |20l3) showed, based on 
modeling the evolution of the luminosity and sizes of high- 
redshift Lyman break galaxies, that at z > 4 the global gas 
accretion rate exceeds the star formation rate, while below 
that redshift they track each other. 

In summary, the evolution of gas fractions reflects a 
competition between cosmic inflow and gas consumption 
rates. A quick estimate of their scalings shows that cosmic 
inflow abates faster than gas consumption, resulting in drop- 



ping gas fractions with t ime in all models. A nalytic models 
based on this scenario bv lDutton et al.l (|201(]| ') show a similar 
result. Outflows provide higher-order modiflcations to this 
picture, particularly owing to wind recycling and preventive 
feedback effects at late times. These trends can be under- 
stood by considering the dependence of the speciflc SFR 
(i.e. tsp) and the depletion time tdcp on mass and redshift. 
At sufficiently early epochs, infiow will be so rapid that gas 
processing cannot keep up, resulting in a gas accumulation 
phase. 



5 COMPARISONS TO DATA 

While the main purpose of this paper is to understand how 
galactic inflows and outflows impact the gas and metal con- 
tent of galaxies, it is instructive to see how our suite of 
models compares to key observations of these quantities out 
to high redshifts. We have already seen that different out- 
flow models generate signiflcantly different predictions for 
gas and metal content. Here we compare to a recent sam- 
ple of forefront observations to see how they constrain our 
outflow models. 

Figure [7] shows a comparison of the mass-metallicity 
relation (top panels) and the mass-gas fraction relation 
(bottom panels) to observations (in black and cyan) at 
z = 0, 2, 3. The z — m etallicity observations are taken 
from pTremonti et al. I (I2OO4I). and the g as fra ction data from 
a compilation bv ^Peepl es fc Shankail (|201C1). At z = 2, w e 
show metallicites and gas fractions from lErb et"al] (|2006al '): 
the gas fractions here are inferred from the star formation 
rate surface density and assuming the Kennicutt-Schmidt 
Law. We also show (in cya n) "direct" gas fraction measures 
from lTacconi et al.l (|20ld ) using CO measurements plus an 
assumed c onversion of CO t o H2 - At z = 3, we compare to 
data from iMannucci eTed} (|2010l ). which is from a sample 
of L yman break ga l a^xies at z ^ 3 — 4 (black points), and 
from iRichard et al.l l|201ll) from a sample of lensed galaxies 
at z = 2.5 — 3.1 (cyan points, metallicities only). 

Looking at the z — MZR, the only model that is 
clearly discrepant is the one with no winds. Galaxies are 
(not surprisingly) over-enriched, because this model greatly 
overproduces the stellar content of galaxies (see e.g. Pa- 
per I). Other works have argued that the z = MZR 
can be r eproduced by varying the ISM star formation ef- 
ficiency (|Tassis. Kravtsov. fc Gnedin 20081: Mouhcine et al.l 
l2008h . or the IMF (jKoppen. Weidner. fc Kroupall2007l ). or 
that hierar chical assembly naturally leads to this shape of 
the MZR (|De Rossi. Tissera. fc Scannapiecol bOOTl ). In our 
simulatio n s, the IMF is always assumed to follow that of 
IChabriej l|2003l ). and the gas depletion timescale tdcp is 
essentially the same in all our models including no-wind. 
Hence in our models, these factors are not responsible for 
the differences in the MZR. In our equilibrium model (eq.[7|, 
the no-wind model overproduces metals because it does not 
eject ISM material to suppress star formation. This in some 
sense represents a "cosmological" star formation efficiency, 
i.e. the amount of cosmic i nfall that is convert ed into stars. 
In our models (Uke those of lBrooks et al.ll2007l '). the MZR is 
governed by the cosmological star formation efficiency rather 
than the ISM one. 

The wind models are generally all in the ballpark of 
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Figure 7. A comparison to observations for the mass-metallicity (top panels) and mass-gas fraction (bottom panels) relation in our 
four wind models at 2 = 0,2,3 (left to right panels). Lines with error bars show the running med ian and la variance within mass bins 
from our r48n384 runs. Data point s show the MZR at z = from SPSS llTremonti et aLlbooi. solid line with 16% - 84% en c losing 
errorbars), at z = 2 from lErb et al.l (jiooii) and at 2 = 3 from iMannucci et al .1 l|20ld. z ^3 — 4; black points) and iRichard et al.l l|201ll . 
z ~ 2.5 — 3.1; cyan points). The momentum-driven wind scalings case comes closest to matching at both z = and z = 2 of the models 
here. For gas fra ctions, z = data is shown from a compilation by Peoples & Shankar (2010|), at z = 2, indirect ga s fractions from 
lErb et ahl l l2006bl . black points) and direct (CO-based) gas fractions from lTacconi et al.l j201(]| . cyan"), and at z = 3 from IMannucci et al.l 
hold ). The momentum-driven wind scalings provide the overall best match, but are still discrepant at < 10^'^ Mq; note that current 
samples measuring gas fractions generally do not include gas-poor galaxies that are often (in our models) satellites. 



the z — SDSS data, but they show non-trivial discrepan- 
cies. At the massive end, all models overproduce the metal- 
hcities in M« > M* galaxies. Recall that we have arbi- 
trarily scaled all our metallicities to match observations at 
M« ~ 10^" Mq, so we cannot improve overall agreement by 
rescaling our yields or metallicity indicators further. One 
possible reason is that there is another feedback mechanism 
that kicks in at M* > AI* . As we saw in Paper I, this is 
already requ ired for supp ressing star formation in massive 
systems (e.g. iGabor et al.|[201L ): it would not be surprising 



if it also suppressed metallicities. We note that the observed 
SDSS sample is solely emission-line galaxies, so does not 
incl ude "red and dead " systems; nevertheless, as shown in 
e.g. ISalim et al.l (|2007l ) and Paper I, even star-forming sys- 
tems seem to have suppressed star formation at the massive 
end compared with a simple extrapolation from lower-mass 
systems. Another possibility is that wind recycling is bring- 
ing in highly enriched material that is elevating the galaxy 
metallicities. Hence this discrepancy may suggest that the 
impact of wind recycling is overestimated in these Simula- 
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tions. In any case, quantitatively matching both the mass 
function and the MZR plateau may provide interesting con- 
straints on such feedback mechanisms. 

At the low-mass end, all wind models remain within the 
observed la variance among SDSS galaxies down to the low- 
est masses probed here (M* ~ 10^ Mq). However, the equi- 
librium model predicts that the simulations with constant 77 
should have a flat MZR at low masses (eq.O. This flattening 
trend is clearly visible in the cw run, and less so in the sw 
case since it kicks in at lower masses. Unfortunately we do 
not yet have the computing power to resolve smaller systems 
within a representative z — cosmological volume to con- 
firm this prediction directly, but higher-resolution r uns done 
to z = 2 do show this trend (|Finlator fc DavelbOOsI '). Obser- 
vationally, iLee et aT I l|2006l ) has found that the MZR contin- 
ues with a similar slope down to M« ~ W^Mq. This would 
clearly be discrepant with expectations for the constant-77 
cases, and would favor models where there is progressively 
higher at lower masses as in momentum-driven wind scal- 
ings. 

The mass-dependent MZR evolutionary rates among 
wind models offer an opportunity to constrain such mod- 
els by comparisons with higher-z data. The no wind model 
continues to over-enrich galaxies by ~ x2 — 3, but the slope 
is now in better agreement. It would be suprising if uncer- 
tianties in yields, metallicity measures, etc. could conspire 
to produce a x2 — 3 change in metallicities, but we cannot 
rule out that possibility. But we disfavor this idea because 
the no-wind model also grossly overproduces the amount of 
stellar mass at z = 2 (as at all redshifts; see e.g. Paper I), 
so it would be highly surprising if it produced the proper 
metal content in galaxies. 

In terms of wind models, at z — the sw model was 
a reasonable fit, but at « = 2 it substantially over-enriches 
galaxies. Since there are few red and dead galaxies at this 
epoch, over-enrichment at the massive end is likely to be 
a serious defect that will not be alleviated by quenching 
(e.g. AGN) feedback. The cw mode l also produces the wron g 
shape for the MZR, as discussed in lFinlator fc Pavel (|2008l ). 
Recall that the MZR shape is probably the most robust pre- 
dictions of our simulations, so again this is a serious defect. 
Meanwhile, the momentum-driven scalings case is a reason- 
able match in terms of shape and amplitude (modulo that 
we have scaled all our metallicities in all wind models by a 
factor of 0.8). 

A t z = 3, the metallicity data from iMannucci et al.l 
lies well below all model predictions, while the data 
from iRichard et all (120111 ) is in good agreement. Since the 
former sample averages a slightly higher redshift, this might 
imply very rapid evolution in the MZR during the short 
time from z ~ 3.5 — > 2.5, but this is very difficult to recon- 
cile with the observed (and predicted) slow evolution from 
z ~ 2.5 0. Models can usually match data at two of the 
three epochs, but no model that we are aware of can match 
all three epochs when considering the lMannucci et al.l (|2010l ) 
data. One possibility is that there are obs ervational selec- 
tion e ffects coming into play, because the IMannucci et al.l 
l|2010f ) z > 3 sample co nsists of rest-u l travio let (UV) se- 
lected galaxies, while the lRichard et all |20l3) sample con- 
sists of lensed galaxies. Rest-UV selection will tend to pick 
out high-SFR systems, which in turn have low metallicities 
for their M« (see Figure [1] with a more detailed discussion 



in It remains to be seen if selection effects can quan- 
titatively reconcile the models and data. Another possibil- 
ity is that there are differences in the poorly- understood 
calib rations of metallicity indicators (e.g. iKewlev fc EllisonI 
l2008h . although both works claim to have carefully consid- 
ered this. Investigating the origin of the differences in the 
observations is beyond the scope of this work, so we merely 
take these as reflective of current systematic uncertainties in 
high-z MZR measures. Upcoming surveys such as the Spitzer 
Extragalactic Deep Survey (SEDS) and the Cosmic Assem- 
bly Ne ar-infrared Deep Extragalactic Legacy Surv ey (CAN- 
DELS; [GroginelaD[20i3;|Koekemo^^ wiU yield 
many stellar mass-selected galaxies at these redshifts, offer- 
ing opportunities through follow-up spectroscopy to measure 
metallicities in large homogeneous samples. 

Turning to gas fractions, at z = all models produce 
the observed trend of lower gas fractions in more massive 
systems. As discussed in i)4.3l the trend in massive (i.e. low- 
/gas) galaxies is set by the depletion time times the specific 
SFR; these both drop with mass, albeit mildly. The gas frac- 
tions in the no- wind case follow the observed shape but are 
too low. One might try to reconcile this by noting that the 
model definition of gas fraction as all star-forming gas (i.e. 
gas above n = 0.13 cm~^) may not be directly compara- 
ble to observational measures of /gas- But at least at face 
value, the no-wind model appears to consume too much of 
its gas into stars by z = 0, leaving galaxies too gas-poor and 
consistent with its overproduction of stars and metals. 

At z = 0, the wind models are broadly in the range 
of observations at high masses, but in all cases they 
show a turnover in gas fractions at low masses that is 
clearly in disagreement with observations. Hence something 
in the current wind simulations is either removing too 
much gas from dwarfs, not supplying enough gas to them, 
and /or consuming their gas t oo qu ickly. Simulations by 
iKobayashi. Springel. fc Whitd (|2007l ') and iMouhcine et all 
1 20081 ) likewise find that dwarf galaxies generally tend to 
be too old with too little present-day star formation, so 
this seems to be a rather generic problem in galaxy for- 
mation simulations within a hierarchical paradigm. In Pa- 
per I we suggested that one explanation may be that the 
star for mation law is different in th ese systems (as argued 
by e.g. iRobertson fc Kravtsovl [20081 ) . particularly at low-z 
when they are typically fairly low surface brightness objects, 
or that the c onversion of neutral into molecula r hydrogen is 
less efficient l|Krumholz. Lerov. fc McKe3l201ll ) . It is curious 
that the simulation without winds at least generally shows 
the correct trend to the lowest masses, and hence another 
possibility is that preventive feedback effects owing to winds 
are incorrectly modeled in these outflow simulations. 

There are other possibilities for the low-mass /gas dis- 
crepancy. For instance, there are more quenched satellites 
at low masses (Figure 4 of Paper I) which have been de- 
pleted of gas. However, as we will show in tj6.3l even central 
dwarf galaxies show a turnover in /gas. Finally, there may 
be observational selection effects in the data as gas mass 
observations (particularly of small systems) tend to focus 
on gas-selected (e.g. H 1) galaxies. For instance, to reconcile 
the vzw model at the lowest mass bin (M« ~ IO^A/q) re- 
quires that observed galaxies are typically 2a outliers in gas 
fractions, which is not impossible. Upcomi ng surveys such 
as the Galex Arecibo SDSS Survey (GASS; ICatinella et all 
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that measure the gas content in an unbiased sample 
will be helpful, but are not targetted to the lowest mass sys- 
tems where the largest discrepancies arise. The gas content 
of low-mass galaxies is therefore an important barometer for 
models, and further data and modeling are needed to shed 
light on whether this is indeed a serious discrepancy. 

At higher redshifts, gas fraction measurements become 
quite uncertain. At 2; ~ 2, all models are in fair agreement 
with observations, although the data overall tend to show 
somewhat higher gas fractions. Selections effects again may 
play a role, since these galaxies tend to be selected as having 
either high gas content (so that direct measures are feasible) 
or high star formation rate (which implies hig h gas con- 
tent; F igure [3l). Indirect gas fractio ns as used bv lErb et al.l 
l|2006al ) andE annucci et all (|2010h have additional system- 
atic uncertainties associated with the validity of using the 
Kennicutt-Schmidt law to infer gas content, while direct gas 
measures must use an uncertain conversion between CO and 
H2 mass. Combined with the broadly similar /gas predic- 
tions among the models, this precludes any meaningful con- 
straints as of yet from high-z gas fractions. The overall rate 
at which gas fractions go down with time at a given M* is 
consistent with observations in all models (modulo at low 
masses at 2 = 0), suggesting that the basic physics govern- 
ing gas fraction evolution has more to do with cosmological 
infall that is same among all models as argued in i )4.3l as 
opposed to feedback mechanisms. 

In summary, all simulations with outflows qualitatively 
reproduce the observed evolution of metallicity and gas frac- 
tion in galaxies. At a given A4"«, all models produce rising 
metallicities and falling gas fractions with time. Quantita- 
tively, the momentum-driven wind scalings model appears to 
be the best overall fit to the ensemble of observations. Nev- 
ertheless, significant discrepancies remain even in this case, 
particularly at the lowest and highest masses, which may be 
owed to systematic uncertainties in metal and gas measures, 
observational selection effects, and model deficiencies. 



6 SCATTER & ITS SECOND PARAMETER 
DEPENDENCES 

The small scatter in the mass-metallicity relation is a strong 
constraint on galactic chemical enrichment scenarios. While 
models have historically focused on reproducing the shape 
of the MZR, the reason for the small scatter, approximately 
0.1 dex v arying only mild ly over more than 5 dex in stel- 
lar mass l|Lee et al.H2"006l ). remains unclear. Models invok- 
ing starburst-induced galactic outfiows would predict larger 
scatter in lower-mass galaxies owing to the more stochastic 
nature of ejection from these systems. 

As is evident from Figure [21 the situation in simulations 
with outfiows is more complex; for instance, our momentum- 
driven wind scalings model produces an MZR scatter at 
2 = that is in good agreement with data, and other mod- 
els less so. Here we examine what governs the scatter in the 
MZR and M* — /gas relations within the context of our equi- 
librium model, particularly focusing on second-parameter 
correlations of metallicity and gas fraction with star forma- 
tion rate and environment. 



6.1 Star formation rate 

The mass-metallicity relation is observed to have second- 
parameter dependences (i.e. correlated scatter) with other 
galaxy properties. A particularly strong one that has been 
explored recently is the sec ond-parameter depen dence of 
MZR on star formation rate. lEllison et al.l (|2008al ) noticed 
that galaxies with high er star formation rates tend to lie 
below the mean MZR. Im annucci et al.l ()2010h developed 
this idea further, and noted that a specific combination of 
Mt and SFR, which they called the fundamental metallic- 
ity relation (FMR), led to a significantly smaller scatter. 
iLara-Lopez et all (|2010l ) independently found a similar re- 
lation. Furthermore, this relation appears to be invariant in 
redshift out to 2 ~ 2.5. Hence the relation between mass, 
metallicity, and SFR provides an even more stringent test 
of model predictions. 

Figure [1] shows galaxies color-coded by star formation 
rate into high (blue), medium (green), and low (red) SFR 
relative to the mean within each stellar mass bin. It is a 
generic result from all simulations that high SFR galaxies 
tend to lie below the mean MZR, and low SFR systems lie 
above. Hence the observed second-parameter trend of the 
MZR on SFR is naturally reproduced in all simulations. 
This suggests a fundamental process at work, independent 
of feedback, that yields this trend. We will now argue that it 
is a natural outcome of the equilibrium model for the MZR 
discussed in i]3.1l 

According to the equilibrium model, a galaxy of a given 
mass prefers to be at an equilibrium metallicity that is estab- 
lished by its current mass loading factor (Equation[7]), which 
sets the amount of enrichment relative to fresh infall. When 
perturbed from its equilibrium metallicity owing to, say, a 
merger, this increases the galaxy mass while lowering the 
metallicity (since smaller galaxies have lower metallicities). 
This moves the system below the mean MZR. Concurrently, 
it results in a more gas-rich galaxy (i.e. a high outlier in 
the MGR), which in turn stimulates star formation. This 
correlated behavior is the origin of the trend that high-SFR 
galaxies lie at lower metallicity. 

Conversely, a galaxy that suffers a lull in accretion will 
consume its gas into stars and form more metals, enriching 
itself along the locus of Z cx M, (when the metallicity is 
significantly below the true yield). This is steeper than the 
MZR slope, and thereby the galaxy moves above the mean 
MZR. Once accretion restarts, the metals become diluted, 
and the galaxy is able to return to the equilibrium relation. 
In some cases, the galaxy is a small satellite dwarf that is 
being environmentally quenched, and it will never restart ac- 
cretion. In the case the sys tem can end up quite far ab o ve the 
MZR, as observations by IPeeples. Pogge. fc Stanekl (120081 ) 
illustrate. 

From the discussion above it is clear that this scenario 
will yield second parameter trends of gas fraction with SFR 
as well. As seen in Figure [S] galaxies with higher SFRs at 
a given mass tend to have high gas fractions. Such galax- 
ies are undergoing an enhanced rate of star formation rel- 
ative to a typical galaxy along the M, — SFR relation, and 
hence will soon consume the gas and return to the mean re- 
lation. In this way, galaxies wobble around the mean MZR, 
M, — /gas, and main sequence relations owing to fiuctuations 
in accretion. Deviations from equilibrium tend to return a 
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Figure 8. The fundamental metallicity relation, Zgas vs. M»SFR"'^^, in our r48n384 runs with four wind models at 2: = 0, 1, 2, 3. Lines 
with error bars show the running median. Small panels below each main panel s how la variance. The v ariance is smaller for the FMR 
as compared to the MZR (Figure [TJ. For comparison, observations of the FMR bv lMannucci et al.l ll2010h found a variance of cr = 0.053. 



galaxy to equilibrium, which is why this scenario is dubbed 
the "equilibrium model." While outflows govern the over- 
all shape and amplitude of these relation, the qualitative 
second-parameter trend with SFR is not a consequence of 
outflows, but rather of equilibrium. This is why the no- wind 
model shows the same second-parameter trend as the wind 
models. 

Outflows do, however, play a significant role in quan- 
titatively establishing the amount of scatter of the MZR 



and MGR. A consequence of t he eq uilibrium model, as em- 
phasized in iFinlator fc Davg (|2008l ) , is that the amount of 
scatter reflects how fast a galaxy can return to equilibrium 
after it suffers a perturbative event. To do so, there must 
be sufficient infall to re-equilibr ate the galaxy. This ca n be 
quantified by the dilution time l|Finlator fc Davil2008l ). 

^^ = {1 + V)-'^ = (1 + V)-'U.„ (10) 

IVJ. il\ IVl if: 

where we have used Equation [6] and the definition of tdcp 
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(eq. El). So long as t^u < tvir, where tvir is the dynamical 
time at the halo virial radius, the scatter about the equi- 
librium relation will be small (|FinIator &: Pavel I2OO8I ). As 
shown in Figure |4l t^cp becomes larger at smaller masses. 
If rj ^constant as in our cw and sw simulations, the scat- 
ter will eventually rise significantly at low masses when the 
dilution time becomes longer than the virial time. This can 
be seen at the lowest masses in the right panels (cw and 
sw) of Figure [2] Conversely, in our momentum-driven wind 
scalings case, the trend of rj with M* mitigates the increase 
in idep to low masses, and keeps the scatter small at smaller 
masses. Specifically, the vzw model (similar to all models) 
yields tdep oc M^"'^, which at 77 2> 1 (small masses) is al- 

— 1/3 

most exactly cancelled by 77 oc M, in this model. So one 
still expects the dilution time, and hence the scatter , to be 
constant to low masses. Observations by iLee et all (I2OO6I ) 
indicate a fairly constant scatter to low masses, though se- 
lection effects may be artificially lowering this (H. J. Zahid, 
priv. comm.). In any case, the scatter in these relations pro- 
vides an independent avenue to constrain ri{M,). 

We now examine more quantitatively the second- 
parameter dependence of the MZR on SFR. Figure |5] shows 
the FMR for our simulated galaxies, namely gas-phase oxy - 
gen abundance versus M*SFR"°-^^ (|Mannucci et al.ll20ia ). 
with the scatter indicated in the smaller panels below each 
main panel. As can be anticipated from Figure[T] the scatter 
is lowered using this combination for the a;-axis. It is not as 
low as observed (a — 0.053 dex), but a somewhat different 
combination of Af, and SFR can lower the model scatter 
further. 

Furthermore, the simulated FMR shows significantly 
less evolution from z = 3 — >■ than the MZR (note that 
the y-a xis scale is signif i cantly smaller than in Figure [2]). At 
z > 3, iMannucci et al.l (120101 ) observes the FMR to evolve 
strongly, which is a consequence of the low observed metal- 
licities in z > 3 _galaxies (see Figure [7]), but observations by 
I Richard et al.l (|201lh indicate no evolution out to z ~ 3. The 
lack of evolution has been taken to indicate that the FMR 
has a special significance that transcends cosmic epoch. In 
our models, however, the lack of (or slow) evolution in the 
FMR is mostly a coincidence; it just so happens that the in- 
crease in metallicity from high-z to low- 2; is balanced by the 
evolution in typical star formation rates for that particular 
combination of parameters. The fundamental principle that 
drives the FMR at any given epoch, namely the tendency 
for galaxies to be drawn towards an equilibrium MZR, has 
little to do with the overall evolution of the SFR at a given 
mass, which is set by cosmic inflow (e.g. Paper I). 

Other observ ations have noted lo wer metallicities in 
merging systems. lEllison et al.l l|2008bl ) showed that close 
pairs tend to have a lower metallicity for their mass. It is 
possible that this trend is driven by the increase in SFR 
in these systems as driven by the interaction, although the 
lack of resolutio n in our models precludes us from examin- 
ing this directly. iPeeples. Pogge. fc Stanekl l|2009l 'l similarly 
found that strong outliers below the MZR tend to be inter- 
acting galaxies and are often quite massive. Once again this 
could be related to their SFR, although massive interacting 
galaxies often have significant AGN activity that produces 
a harder radiation field within their ISM, which when using 
abundance ratios to measure abundances can mimic a lower 
metallicity (C. Tremonti, priv. comm.). While these trends 



are interesting, it is unclear whether they are distinct from 
the overall trend of lower metallicities in high er SFR galax- 
ics (at a given M,). Finally, we mention that Elliso n et al.l 
(2008a) found that systems with large half-light radii, like 
with high SFR, also lie below the mean MZR. Unfortunately, 
our simulations lack sufficient resolution to robustly model 
galaxy sizes, so we cannot directly examine this second- 
parameter dependence. 

The scatter in /gas should, according to this scenario, 
follow the same trend as for the metallicity. Figure [6] s hows 
that in the compilation of IPeeples Sz Shankaij ([2013), the 
scatter is lower at small masses. This may be partly a se- 
lection effect, since particularly at small masses H I or CO 
selected samples will pick out the most gas-rich systems, 
artificially loweri ng the scatter. I n the stellar mass-selected 
GASS sample of ICatinella et al.l (|2010l ). there is no obvious 
evidence for a change in scatter in /gas versus M, among 
the galaxies that contain gas (excluding gas-poor passive 
systems), mirroring the roughly constant MZR scatter. In 
the context of the equilibrium model, this again argues for 
a mass loading factor that increases to lower masses. 

In summary, galaxies in our simulations with high star 
formation rates at a given mass also have lower metallicities 
and higher gas fractions. This second-parameter dependence 
of the MZR and MGR is a natural and straightforwardly 
understood consequence of the equilibrium model. This de- 
pendence is independent of outflows, and arises purely as a 
consequence of equilibrium. Quantitatively, the dependence 
of scatter on M, is a function of rj and tdep, and mod- 
els that have larger mass loading factors at lower masses 
better match observations of constant or mildly increasing 
scatter in the MZR and MGR to the lowest masses. The 
FMR provides an interesting tool to quantitatively examine 
the second-parameter trend with SFR, and simulations yield 
broadly similar trends to those observed for the FMR. 

6.2 Environment 

Another second-param eter depe ndence of the MZR was 
noted bv i Cooper et al.1 (|200a ') and iEUison et all (|2009l '). who 
showed that galaxies within dense environments such as 
groups and clusters tend to have higher metallicities. Galax- 
ies outside such environments, in contrast, tend to have no 
obvious dependence on local galaxy density. This trend is 
over and above any trend associated with star formation. 
Our 48/!.~^Mpc volume has some galaxy groups up to virial 
masses of ~ 10" Mq, but not enough to compare directly to 
observations of clusters. However, we can examine the trend 
in metallicity with environment as measured by local galaxy 
density. 

Figure[H]shows the z — MZR for our four wind models, 
where we have subdivided galaxies by local galaxy density as 
measured in a l/!,~^Mpc tophat sphere. Galaxies at densities 
> 0.5a" above the mean are shown in red, < 0.5a below the 
mean in blue, and those in between in green. A running 
median for the MZR, with la variance, is shown for each 
sub-population. 

The momentum-driven scalings and no-wind cases dis- 
play the observed trend: Galaxies in high-density regions 
lie above the mean MZR by ~ 0.05 dex, while galaxies in 
medium and low density regions show no discernible differ- 
ence. The constant-77 models show signiflcantly less depen- 
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dence on environment. The difFerences disappear at the most 
massive end in aU models. 

One reason for this dependence may be that denser 
environments have more e nriched intergalactic gas (e.g. 
lOppenheimer fc Pavel [200^ ). so accretion onto galaxies in 
those environs is likely to boost the metallicity over galaxies 
in a less dense region. We believe this is indeed the trend 
responsible, but how this operates is subtle, and gives rise 
to distinct features among the wind models. 

The model trends are best understood if the metallic- 
ity of infalling gas is governed by wind recycling. Recall 
that the IGM is almost e ntirely enriched by winds (e.g. 
lOppenheimer fc Davell2006l ). and so any metals falling back 
into galaxies constitutes wind recycling. At high masses, 
wind recycling is so effec tive that all galaxies re-acc rete their 
ejected material quickly (|Oppenheimer et al.ll2010l ). which is 
like having no winds at all, so the metallicity approaches the 
overall yield regardless of environment (modulo an increase 



due to enriched infall; see JT]). At sufficiently low masses, 
all ejected material escapes, and hence the infalling mate- 
rial is mostly primordial regardless of environment. But in 
the intermediate regime, environment plays a critical role 
in slowing winds (|Oppenheimer fc Davell2CI0g '). in the sense 
that denser regions slow winds more and cause faster recy- 
cling. 

As shown in 'Oppenheimer et al.l (|2010t ). in the case 
of momentum-driven scalings this intermediate regime oc- 
curs over a protracted range in M* since oc v^sc', this 
protracted mass range is reflected in this model's depen- 
dence of MZR on environment. In the constant-v^, cases 
(cw and sw), there is only a small range of masses between 
the fully-escaping and fully-recaptured regimes (around the 
mass where « Vcbc), which means that the environmental 
dependence is only seen over ~ 0.5 dex in mass. Hence in 
these models, the high-density MZR actually shows a peak 
in metallicity at the mass where recycling is most effective. 
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It occurs at higher in the cw case relative to sw since its 
higher wind speed allows escape up to larger galaxies. 

Empirically in our simulations, the dependence of 
metallicity on environment seems to be only effective 
when the environment becomes quite dense; medium and 
low density regions show no difference. This is likely be- 
caus e only these dense environments have hot gaseous ha- 
los H eres et al ] l2009a') that can significantly slow winds. 
The no-wind case also shows a dependence on environment 
that is obviously not driven by wind recycling, but may be 
driven by tidal stripping (and subsequent enriched infall) 
which is more effective in dense regions; indeed, in 33 we 
show that the no- wind case at z = has non-negligibly en- 
riched infall. We leave a more detailed study of these effects 
for the future. Here we simply suggest that the environmen- 
tal dependence of wind recycling governs how local galaxy 
density impacts the MZR, and our momentum-driven wind 
scalings model generally reproduces the observed trend. If 
this is true, studying the environmental dependence of the 
MZR offers a unique probe into cycle of baryons in and out 
of galaxies. 

6.3 Satellites 

Satellite galaxies are seen to have higher metallic ities at a 
given stellar mass (|Peeples. Pogge. fc Stanekl2009l ). The ori- 
gin of this trend is qualitatively understood within the equi- 
librium model: Satellites tend to live in dense regions where 
the infalling gas is more enriched, and furthermore they are 
not straightforwardly fed by cold streams since they do not 
lie at the bottom of the halo's potential well and hence evolve 
upwards off the MZR. In this section we quantitatively assess 
the differences in satellite versus central populations for the 
various second-parameter trends we have examined above. 

Figure [10] shows the median mass-metallicity relations 
for a variety of second-parameter dependences. In order to 
more easily see the dependences, we have subtracted out the 
overall MZR from each subsample's MZR. The plot shows 
A[O/H]0 dependences on three variables: Star formation 
rate, environment, and satellite vs. central galaxies. For the 
first two quantities, we subdivide the overall sample into 
"high" and "low" , which simply means above and below the 
median within each mass bin. In the left panels, we show 
the median (differenced) MZR for galaxies with high-SFR as 
blue, and low-SFR as red. In the right panels, we analogously 
show the MZR for high-density environment galaxies in red, 
and low in blue. Top panels show the vzw simulation, and 
bottom panels show the no-wind run. 

We further explore the dependence of SFR and envi- 
ronment within satellite (dotted lines) and central (dashed) 
galaxy samples. This is done for the full sample (black), 
high-SFR/low-density (blue) and low-SFR/high density 
(red). The solid line at represents the original MZR of 
all galaxies. In all, this figure shows how the dependences 
on SFR and environment interplay with the distinction be- 
tween central and satellite galaxies. 

Let us begin examining Figure [TO] by considering the 
blue and red solid curves, i.e. the MZR subdivided by SFR 
and environment. These trends have been noted earlier, but 
here the differences are more visible since we have sub- 
tracted off the overall trend. Comparing the solid red and 
black curves, we see that, as before, low-SFR and high- 



environment galaxies lie above the global MZR, at least for 
galaxies with M* < 10^^ Mq. The trend with environment 
is generally stronger than the trend with SFR, which may 
be surprising given that this second-parameter dependence 
has received less attention in the literature. In the context 
of the equilibrium model, this indicates that denser envi- 
ronments result in significant suppression and enrichment 
of inflow, such that galaxies within dense regions both have 
lower star formation rates (causing higher metallicities) and 
are accreting higher metallicity gas. 

The source of these trends is further clarified when 
examining the satellite galaxy population. Overall, satellite 
galaxies (dotted curves) show elevated metallicities at a 
given mass compared to central galaxies (dashed curves), 
mimicking observed trends (jPeeples. Pogge. fc Stane^ 
[200^. Since satellites have typically lower SFRs (Pa- 
per I) and also tend to live in denser environments, both 
second-parameter dependences discussed previously could 
be contributing the satellites' higher MZR. A particularly 
striking result is that the second-parameter dependence on 
SFR is fairly small in the central galaxies, and is dominated 
by the satellite systems. That is, there is an enormous 
difference between the metallicities of high-SFR satellites 
(dotted blue) versus low-SFR satellites (dotted red) over 
most of the sub-M* mass range. The trend is very strong in 
the vzw run, but also present in the nw run. The detailed 
physical origin of this we leave for future work, but here 
we speculate that galaxies first entering another halo will 
have enhanced SFR due to interactions and therefore will 
have lower metallicities, but will then get quenched and 
quickly build up metallicity to go above the mean MZR. 
In any case, it is evident that in these simulations, the 
second-parameter dependence of the MZR on SFR is driven 
more by the satellite galaxies, while centrals show only a 
modest such dependence. 

Looking at satellites vs. centrals subdivided by envi- 
ronment, it is not as clear here what is driving the overall 
trend. Both centrals and satellites show higher metallicities 
in dense regions. In the vzw model, the trend is somewhat 
stronger for satellites, but not as dramatically as in the case 
of SFR. Satellites are affected by both diminished (from 
strangulation) and enriched inflow, while central galaxies 
should not have inflow enriched, only diminished since they 
should still accrete gas normally at the bottom of the halo's 
potential well regardless of environment. That the strength 
of the affect is only slightly stronger for satellites indicates 
that the majority of the change to the MZR in dense regions 
arises because inflow is more enriched owing to residing in a 
denser region. It is worth noting that these changes in metal- 
licity are, in an absolute sense, not large: they are typically 
below 0.1 dex. There fore even a modest metall icity-density 
gradient in the IGM (|Qppenheimer et al.ll201ll ) could pref- 
erentially cause galaxies in denser environments to be over- 
enriched by this amount. 

In summary, both environment and star formation rate 
play an important role in driving second-parameter trends 
in the MZR. The majority of this trend is driven by satellite 
galaxies, as central galaxies are less affected by these second 
parameter trends. This shows, as expected, that satellites 
have greater fluctuations in their accretion rates that drive 
star formation, and are more impacted by environmental 
effects. The overall trends are consistent with expectations 
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Figure 10. Second-parameter dependences of the MZR subdivided by satellite vs. central galaxies at z = in our simulations. Each 
curve shows the difference between the MZR for a particular subsample of galaxies; the solid black curve at zero shows the overall 
MZR for that model. Top panels show the momentum-driven wind scalings case, bottom panels the no-wind simulation. Left panels 
show galaxies subdivided by star formation rate (at a given mass) above and below the mean. Right panels show galaxies analogously 
subdivided by environment. The blue curves show galaxies with above-median SFRs (in left panels) and below-median environments 
(right panels), red curves show below-median SFRs and above-median environments. Dashed and dotted black curves show the MZR for 
satellites and centrals, respectively. These populations are further subdivided into high-SFR/low-environment (dashed and dotted blue 
curves) and low-SFR/high-environment (dashed and dotted red curves). 



from the equilibrium model, being driven by a competition 
between recent accretion and outflows. In the case of satel- 
lites vs. centrals at a given stellar mass, the outflow rates 
are similar in the simulations (being zero in the no-wind 
case), but the accretion rates can vary owing to both gen- 
eral stochastic fluctuations and environmental suppression. 
The typical difference in metallicities between satellites and 
centrals is quite small, typically of order 0.1 dex, so par- 
ticular care is needed to tease out such effects in observed 
samples. 



7 EVOLUTION IN GAS & METAL CONTENT 

We have previously studied the evolution of the overall MZR 
and MGR relations, finding that at a given mass, metallici- 
ties rise slowly and gas fractions fall slowly with time. In this 
section we examine in more detail how particular galaxies 
evolve within these relations, in order to better understand 
the nature of the overall evolution. 

FieurelTTIshows galaxy tracks in MZR space (top panel) 
and MGR space (bottom panel) from = 2^0. We choose 
our vzw model since it does the best job of matching observa- 
tions (with notable exceptions) of the models considered. To 
make these tracks, three representative stellar masses were 
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Figure 11. Evolution from 2: = 2 — > of mean metallicity (top 
panel) and gas fraction (bottom) versus stellar mass for a set 
of galaxies within 3 mass bins in our r48n384vzw run. Cyan and 
magenta points show the overall galaxy population at z = 2, 0, re- 
spectively. Numbers along the tracks indicate the redshift; tracks 
end at 2; = 0. 



selected at « = (specifically Vd^'^-^'^ , lO^^-^^lO^'®'^ Mq), and 
20 galaxies were chosen closest to each mass. The main pro- 
genitors of each galaxy were identified in each output back 
to z = 2, where the main progenitor is the galaxy at an ear- 
Uer epoch hosting the largest fraction of the final galaxy's 
particles. The tracks shown are the mean value of the 20 
progenitors. Numbers along the tracks indicate the redshift, 
with the tracks ending in an arrow at z = 0. 

The most evident trend is that galaxies tend to evolve 



mo stly along the mean MZR and MGR relations, as noted 
by iBrooks et ahl l|2007l ) for the MZR. This directly trans- 
lates into a slow evolution for these relations. The fact that 
galaxies move along these relations has sometimes been for- 
warded as the "cause" for the slow evolution, but this merely 
begs the question, why do galaxies tend to move along these 
relations? 

To answer this, let us first consider the MZR. It is 
straightforward to differentiate Equation [7| to show that if 
rj oc M^^ , then dlogZ/dlogM, — x when 3> 1; this is 
the equilibrium model prediction for the slope of the galaxy 
track in MZR space at low masses, assuming that the con- 
stant of proportionality for r\ for a given galaxy is unevolv- 
in^. Hence if a galaxy obeys Equation [7] at all times, it will 
evolve directly along the MZR relation when 77 ^ 1. For 
example in the vzw case, x — 1/3, which is identical to the 
MZR slope in the low-M* regime. 

However, Figure [11] shows that the evolutionary slope is 
steeper than this: for the lower mass bins, d log Z/d log M, « 
0.6. The more rapid evolution must arise because an assump- 
tion in Equation [7] is violated. In particular, it turns out 
that the infall is not pris tine as assumed in that equation. 
lOppenheimer et al.l l|201ll . Figure 3) shows that the metallic- 
ity just outside star-forming regions (i.e. at uh ~ 0.13 cm"'') 
rises substantially from 2 = 2 — J- 0, and exceeds solar today. 
They argue that this arises owing to the preponderance of 
recycled wind accretion at later epochs, which causes the 
inflow from t he IGM to be increasing ly enriched. 

Following iFinlator fc Pavel l|2008h . we can extend Equa- 
tion [7] to include the effects of enriched infall. If we define 
az as the ratio of infalling gas metallicity (^infaii) to the 
metallicity within the ISM of the galaxy (Zism), then 



y 



1 



1 + rj 1 - az 



(11) 



We can directly measure az in our simulations. For a given 
mass, we take 50 galaxies near that mass and compute the 
mean metallicity within the ISM (i.e. star-forming) gas. We 
then compute the mean metallicity in infalling gas. We de- 
fine infalling gas as all gas within 30 kpc (comoving) that 
is not star-forming and is moving towards the galaxy (i.e. 
v-r < 0). We also tried scaling the infall radius with the 
virial radius at different masses, with only minor differences. 

The evolution of Zism and ^infaii are shown as the 
solid and dotted lines in the top panel of Figure 1121 We 
show the vzw model at two masses, and the no-wind case 
at Mt ~ 10^*^ M0. The metallicity is higher around larger 
galaxies, as expected. The interesting trend is that Zinfaii 
increases faster than Zism. This is quantified in the bottom 
panel where we plot (1 — az)~^ , which is the extra factor in 
Equation [11] accounting for enriched infall. The key point is 
that, from 2 = 2-^0, this term increases by 0.2 — 0.3 dex. 
This is identical to the excess increase in the galaxy metal- 
licity above simply moving up along the MZR (Figure [TT|) . 



This turns out to be basically true, although it is something of 
a coincidence: Equation |3] shows that for a given galaxy mass, a 
drops with time, implying a higher 77 and thus a lower metallicity. 
However, this is countered by the fact that a given galaxy's a 
increases with time as it grows. It happens that the two effects 
mostly cancel out for any given galaxy. 
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Figure 12. Top panel shows the evolution from z = 3 — >■ of 
mctallicity in the ISM Zism (solid lines) and metallicity of in- 
falling gas ^infall (dotted lines). Infalling gas is all gas within 30 
comoving kpc that is not star-forming and is moving towards the 
galaxy. Green and cyan curves show results for the vzw simula- 
tion at stellar masses of 2 X 10^ and 10^" Mq, respectively. The 
blue curves show the results for the no- wind case at 10"'^'^ Mq. 
Bottom panel shows the extra multiplicative term in the equilib- 
rium model MZR (Eguation lllH that accounts for enriched infall, 
namely (1 — az)~^ where az = .^infall/'^lSM. This shows that 
enriched infall is primarily responsible for the upward evolution 
of the MZR in our simulations. 



The perhaps surprising imphcation is that the rising metal- 
licity at a given stellar mass is not the result of galaxies 
processing more gas into stars, but rather the result of an 
increasing metallicity in accreted gas. 

The no- wind case also shows a similar increase in (1 — 
az)~^ of about 0.2 dex from z = 2 — > 0. In this case, this 
arises because tidal interactions distribute metals around 
galaxies that can later fall back in. This has only become 
prominent since z ~ 1. Overall, the enriched infall term is 
much lower than in the vzw case, showing that most of the 
enrichment in the infall is generated by outflows. 

The gas fractions of our selected galaxies also show an 
evolution generally along the relation. The evolution does 
become notably steeper at low redshifts {z < 0.5 — 1) and 
at small masses, which shows that the "turnover" in low- A/* 
gas fractions is a late-time phenomenon. These dwarf galax- 
ies are apparently depleting their gas reservoir too quickly, 
and are prevented from re-acquiring their ejected material 
owing to preventive feedback processes. We leave a fuller ex- 
amination of the interplay between such feedback processes 
and gas content in dwarf galaxies for the future. 

In summary, the slow evolution of the MZR and MGR 
scaling relations procedurally arises from the fact that galax- 
ies tend to evolve mostly along these relations, with only 
mild deviations towards higher metallicity and lower gas 
fractions with time. Within the equilibrium model, the 



higher metallicities arise because gas infall is increasingly en- 
riched to lower redshifts, while the lower gas fractions arise 
because the depletion time becomes smaller compared to 
the star formation timescale. While these trends are quali- 
tatively consistent with the idea that gala^xies obtain a large 
reservoir and slowly consume their gas (while generating 
metals), our simulations suggest that the actual physics is 
much more complex, driven by a balance between inflow and 
outflow processes. 



8 SUMMARY 

In th is paper and Paper I (|Dave. Oppenheimer. fc Finlatoi] 
we have presented a study of how the stellar, gas, and 
metal contents of galaxies are governed by gas inflow and 
outflow processes within an evolving hierarchical Universe. 
In Paper I we investigated how galactic outflows play a key 
role in modulating the stellar growth of galaxies fed pri- 
marily by cold, filamentary accretion from the IGM. In this 
paper, we have shown that such infiow and outflow processes 
concurrently govern the evolution of the metallicity and gas 
fraction within star-forming galaxies. 

The central message of these two papers is that the evo- 
lution of the main constituents of star-forming galaxies can 
be broadly understood within the context of a cycle of in- 
flow and outflow between galaxies and the IGM. An idea 
that features prominently in our models for the evolution 
of the gas and metal content is the notion of equilibrium. 
Galaxies prefer to live on specific equilibrium relations be- 
tween metal, gas, and stellar content, whose forms are set 
by the cosmologically-evolving infiow and outflow rates. The 
inflow rate into the ISM is tied to the accretion rate into ha- 
los, with notable departures at low masses and late epochs 
owing to preventive feedback. Meanwhile the outflow rate 
appears to be most closely tied to stellar mass, since e.g. 
the metallicity of a galaxy is observed to have the tightest 
correlation with its stellar mass as compared to any other 
individual property. Inflow fuels star formation, whereas out- 
flows are the central governing agent that control how much 
of the inflowing material turns into stars. An key corollary 
of this scenario is that stochastic variations in the inflow 
rate tend to drive galaxies back towards the equilibrium re- 
lations, resulting in small scatters in metallicities and gas 
fractions that are correlated with star formation rate. This 
equilibrium paradigm can therefore quantitatively explain 
the origin of the shape, slope, and scatter of the relations 
between gas, metals, and stars, as arising naturally from 
hierarchical galaxy growth modulated by outflows. 

With that framework in mind, we summarize the key 
conclusions of this paper: 

• Galaxy metallicities are set by a balance between in- 
flows that provides (relatively) pristine fuel, and outflows 
that reduce the cosmological star formation efficiency by 
ejecting fuel. This equilibrium can be expressed as a func- 
tion of the outflow's effective mass loading factor 77 (eq. [7|), 
meaning that the stellar mass-metallicity relation (MZR) 
mostly reflects the relationship between 77 and stellar mass. 

• The evolution of the MZR in this scenario is expected 
to occur mostly along the relation, as confirmed by tracking 
simulated galaxies. In detail there is a slow upwards evolu- 
tion in metallicity at a given stellar mass, i.e. the MZR rises 
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with time. We demonstrate that this is quantitatively un- 
derstood as a result of accreted gas becoming more enriched 
with time. 

• Galaxy gas fractions reflect a competition between gas 
accretion, as quantified by the star formation timescale 
(Mt/SFR), and gas consumption, as quantified by the de- 
pletion time (Mgas/SFR) (Equation [9l . The depletion time 
in our models is set primarily by our assumed law for star 
formation (based on Kennicutt-Schmidt), while the star for- 
mation time is governed by cosmic inflow. Since both of these 
timescales are relatively insensitive to outflows, gas fractions 
are (in contrast to metallicities) likewise insensitive to out- 
flows. 

• The stellar mass-gas fraction relation (MGR) drops 
slowly with time in all models. This arises because the gas 
supply rate, driven by cosmic accretion, drops faster than 
the gas consumption rate, which is tied to the galaxy's dy- 
namical time. Constant replenishment of ISM gas is a ubiq- 
uitous feature of these models, as appears to be required 
from observations. Galaxies individually evolve mostly along 
the MGR, but drop particularly at late epochs. 

• Wind reaccretion plays an increasingly important role 
in the evolution of the MZR and MGR at late times, par- 
ticularly at 2 < 1. Enriched inflow directly corresponds to 
material that was ejected at an earlier epoch, and alters the 
MZR (Equation [TT} . In the MGR, all wind models develop 
a turnover at late epochs and low masses, likely reflecting 
the lack of re-accretion of wind material in these systems. 

• The scatter in the MZR and MGR reflects how fast 
a galaxy can return to equilibrium given a fluctuation in 
the accretion rate, as quantifled by the dilution time given 
by tdep/(l + v)- The scaling of the MZR and MGR scatter 
with mass therefore provides an independent constraint on 
?7(M.). 

• Departures from equilibrium naturally correlate with 
star formation in all models, even without winds, as it is a 
consequence of equilibrium rather than feedback. Galaxies 
with high SFR for their M« are predicted to have low metal- 
licity and high /gas, consistent with observations. Galaxies 
in denser regions and satellites are also predicted to have 
higher metallicities as observed, since these systems are ob- 
taining more enriched inflow and/or their inflow has been 
curtailed because they are not residing at the center of the 
halo. 

• Comparing to observations of the z = MZR, the equi- 
librium model with momentum-driven wind scalings pre- 

1/3 

diets an unbroken power law of Z oc M* with small scatter 
as broadly observed, while the constant-?7 models predict a 
flattening of the MZR at low masses with a strongly increas- 
ing scatter. All models qualitatively reproduce the observed 
slow evolution of the MZR upwards from z = 3 — > 0, but 
the momentum-driven wind scalings model comes closest to 
matching data at both z = 2 and z = 0. 

• Gas fractions at z = in all wind models predict a 
falling /gas with M* at high masses, and a turnover to lower 
/gas at the smallest masses. The former broadly agrees with 
data, while the latter is in clear disagreement with data. This 
turnover is related to a downturn in speciflc SFR and upturn 
in age in small systems (Paper I), and likewise may indicate 
that star formation in dwarfs must be delayed on cosmic 
timescales, perhaps owing to a different star formation law 
or less efficient conversion of HI to H2 in these systems. 



Considering both metallicities and gas fractions in this 
paper along with stellar masses and star formation rates 
from Paper I, our simulations with momentum-driven wind 
scalings provide the best ensemble match to available ob- 
servations of the models considered here. Its general success 
hinges on having a higher mass loading factor in smaller 
systems, as well as not having a characteristic velocity scale 
picked out by the wind model that would translate into dis- 
tinct features in relations such as the MZR that are clearly 
not observed. In both papers, however, we found that the 
model best matches at mass scales from M* down to mod- 
erately large dwarfs. At larger masses, it is clear that some 
form of quenching feedback is required, which is manifested 
here as overly high metallicities and a dearth of gas-free 
galaxies. At small masses, Paper I showed how small dwarfs 
have too low star formation rates, which is accompanied 
here by too low gas fractions. Nevertheless, this is the first 
cosmological hydrodynamic simulation that provides a rea- 
sonable match to observed present-day galaxies around Ad* 
that contain the majority of cosmic star formation. 

Taking a broader view, the problem of galaxy evolu- 
tion appears to be separable into three phases, predomi- 
nantly divided by halo mass, which could somewhat fanci- 
fully be called the birth phase, the growth phase, and the 
death phase. In the birth phase, halos are small enough that 
photoionization plays a critical role in retarding ac cretion; 
these are galaxies below the so-called filtering mass (i GnedinI 
|2000i). In the growth phase, galaxy growth is regulated by 
baryon cycling, i.e. by smooth, filamentary accretion and 
ubiquitous outflows that circulate mass, energy, and metals 
between galaxies and the IGM. In the death phase, some 
mechanism (often associated with feedback from the central 
black hole) quenches accretion into the ISM, thereby halting 
star formation and creating a passive galaxy. Thus the three 
phases are not only distinct in halo mass but also have three 
separate dominant feedback mechanisms that govern galaxy 
growth, making galaxy evolution in each phase physically 
distinct from the others. 



In this paper and Paper I we have focused on galax- 
ies in the growth phase. In this phase, the dark matter halo 
virial radius as a boundary between the galaxy and the IGM 
is of secondary importance: gas flows in and out unabated 
through the virial radius. Major mergers are a sub-dominant 
fueling mechanism in such systems, and are relatively unim- 
portant in the overall star formation history. Central black 
holes may be growing within these systems, but they play a 
minor role in the overall evolution. In contrast, the formation 
and evolution of passive "death phase" galaxies appears to 
be critically linked to major mergers, black holes with asso- 
ciated fee dback, and the pre sence of a stable hot gaseous 
halo (e.g. iKeres et al.l [2005") that de marcates the virial 
ized region f rom the ambient IG M (e.g. ICroton et al 
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tion between the growth and death phases may be triggered 
by a major merger, such t hat bulge f ormation and black 
hole growth are linked fe.g. iHopkins et al. 200&), although 
it appears that maintaining a red and dead gal axy primarily 
relies on the existence o f a hot gaseous halo (jCroton et al.l 
l2006l : [Cabor et al]l201ll l. While these phases are physically 
distinct, understanding the evolution of the galaxies in the 
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death phase hkely requires weU-known "initial conditions" 
provided by galaxies in the growth phase. 

Many key aspects of the galaxy life cycle are far from 
being fully understood. It will be a transformative achieve- 
ment in the galaxy formation community when even the ba- 
sic framework is in place, and that work can begin towards 
a quantitative rather than qualitative understanding of the 
main physical processes (analogous to the era of precision 
cosmology). In these two papers, we have illustrated how 
simulations can be used to elucidate simple analytic relation- 
ships between physical processes of inflow and outflow and 
the observable properties of galaxies. We have shown that 
cosmic dark matter-driven inflow combined with the mass 
loading factors, wind recycling properties, and preventive ef- 
fects of outflows govern evolution of the basic constituents 
of galaxies. This provides a bridge between detailed studies 
of inflow and outflow, the latter being much more poorly un- 
derstood, and large-scale surveys of galaxy properties across 
cosmic time. In other words, given a detailed model (e.g. 
from individual halo simulations) for outflow mass loading, 
wind recycling, and preventive effects, we can now translate 
that with good fidelity into predictions for the evolutionary 
properties of galaxy populations. Alternatively, it allows ob- 
servations of galaxy populations to be straightforwardly in- 
terpreted as constraints on the detailed properties of inflows 
and outflows. This scenario provides a first step towards 
understanding the much wider range of interesting observ- 
able galaxy properties, such as kinematics, radial gradients, 
morphologies, and environmental dependences. Continuing 
to grow the synergy between multi-scale models and multi- 
wavelength observations of galaxies and their surrounding 
gas is the best way to advance our understanding of the life 
cycle of galaxies across cosmic time. 
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